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Abstract 

The  high  energy  products  of  the  10B(n,a)'Li  neutron  capture  reaction  were  ex¬ 
plored  as  a  means  of  thermal  neutron  detection  in  thin  film  amorphous  boron  carbide 
(a~B5C:Hx)  on  n-type  silicon  hetero-junction  diodes.  The  diodes  were  produced  via 
plasma-enhanced  chemical  vapor  deposition.  A  second,  resistive  type  device  exploit¬ 
ing  the  same  reaction  products  to  detect  neutrons  was  also  tested.  Pulse  height 
spectroscopy  was  conducted  on  the  diodes  while  they  were  subjected  to  a  moderated 
plutonium-beryllium  (PuBe)  source  with  a  thermal  neutron  flux  of  464  [n/cm2-s]. 

Neutron  capture  could  not  be  confirmed  experimentally  in  any  of  the  devices 
tested.  Cadmium  shielding  was  used  to  discriminate  between  neutron- induced  signals 
and  other  interactions.  Both  diodes  failed  to  produce  a  pulse  height  spectrum  with 
identifiable  features,  either  with  or  without  cadmium.  The  activity  observed  from  the 
resistive  device  was  conclusively  determined  to  not  be  a  result  of  a  neutron  interaction. 
The  cadmium  wrap  produced  pulses  with  a  greater  frequency  and  amplitude  than 
previously  observed  in  the  unshielded  spectra,  and  the  observed  signal  was  attributed 
to  photon  down-scattering  in  the  cadmium  shield. 

Device  modeling  was  performed  to  investigate  the  cause  for  the  diode  detector’s 
negative  response.  A  simulated  1.47  MeV  a  particle  produced  transient  current  of 
10"'  A  and  the  induced  charge  was  4.7  times  less  than  intrinsic  noise  charge  of  the 
preamplifier.  Further  investigation  of  the  model  emphasized  the  transient  current 
magnitude  was  most  dependent  on  the  electric  field  strength.  The  possibility  of 
direct-conversion  hetero-junction  neutron  detectors  was  not  experimentally  confirmed 
in  this  study.  Proof-of-principle  simulations  using  only  10  /irn  of  silicon  however  did 
generate  transient  currents  2  factors  larger  than  steady-state  reverse  bias  currents. 
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NEUTRON  DETECTION  USING  AMORPHOUS  BORON-CARBIDE 


HETERO- JUNCTION  DIODES 

I.  Introduction 


1.1  Background 

Through  the  end  of  the  Cold  War  in  1991,  the  stockpiles  of  special  nuclear  material 
(SNM)  were  tightly  controlled  by  the  world’s  most  powerful  states.  Special  nuclear 
material  as  defined  by  Title  I  of  the  Atomic  Energy  Act  of  1954  includes  Pu,  233 U, 
and  uranium  enriched  in  the  isotopes  233U  or  235 U  [1].  The  fall  of  the  Soviet  Union 
precipitated  a  lapse  in  the  security  of  SNM  with  profound  effects,  and  the  threat  of 
proliferation  of  SNM  has  been  a  global  concern  ever  since. 

Several  nations  have  gained  possession  of  SNM  since  1991.  Some  of  these  states, 
especially  the  Republic  of  North  Korea  and  Iran,  have  adversarial  stances  with  the 
United  States.  Nuclear  states  with  unstable  governments  pose  considerable  risks  to 
the  United  States,  and  to  some  degree,  the  rest  of  the  world.  Extremist  organizations 
have  demonstrated  their  desire  to  inflict  harm  on  the  US  and  its  allies;  they  are 
limited  only  by  their  means  to  acquire  the  tools  to  do  so.  The  2001  terrorist  attacks 
against  the  United  States  removed  all  doubt  of  this,  and  for  the  foreseeable  future, 
changed  the  American  sense  of  security  and  made  the  threat  of  SNM  proliferation  a 
relevant  concern. 

Detection  at  standoff  distances  of  a  minimum  of  several  meters  is  an  inherent 
requirement  for  intercepting  any  attempt  to  smuggle  SNM  into  the  US  through  regular 
points  of  entry,  especially  in  shipping  ports  where  large  containers  help  to  cover  up 


1 


any  nuclear  signature.  Detecting  SNM  presents  an  additional  challenge  for  standoff 
detection.  Most  isotopes  of  SNM,  240Pu  being  the  only  exception,  does  not  have 
significant  spontaneous  fission  neutron  emission  rates.  240Pu  has  a  spontaneous  fission 
neutron  emission  rate  of  ~920  [n/g-s]  [2],  Furthermore,  gamma  ray  emission  for 
typical  SNM  isotopes  with  nontrivial  emission  rates  is  mostly  below  200  keV  for 
decays  [3].  The  limited  range  for  gamma-rays  and  neutrons  makes  it  very  easy  to 
shield  highly  enriched  uranium  (HEU)  and  weapons  grade  plutonium  (WGPu). 

Cosmic-ray  induced  spallation  neutrons  further  complicate  matters  since,  when 
generated  at  high  Z/air  interfaces,  can  cause  a  neutron  background  of  '0.077  [n/cm2-s]. 
Typical  background  for  air/sea  interfaces  is  less  (around  0.012  [n/cm2-s]).  A  hypo¬ 
thetical  weapons  mass  of  4  kg  of  94%  239Pu  WGPu  is  detectable  in  an  environment 
where  high  Z/air  interfaces  abound.  This  hypothetical  mass  of  WGPu  has  an  isotropic 
neutron  flux  on  the  order  of  0.39  [n/cm2-s]  at  a  distance  of  3  m  assuming  no  shield¬ 
ing  or  moderation  [4],  A  hypothetical  mass  of  HEU  (25  kg  for  an  implosion-type 
weapon)  on  the  other  hand  emits  neutrons  at  a  rate  four  orders  of  magnitude  less 
than  WGPu,  resulting  in  a  neutron  flux  well  below  background  levels  at  3  m.  If  there 
was  significant  gamma  shielding  but  no  neutron  specific  shielding,  it  is  feasible  that 
WGPu  could  be  detected  from  its  neutron  emission  signature  while  HEU  would  go 
undetected  due  to  its  weak  neutron  emission  activity. 

The  method  of  SNM  detection  described  in  this  thesis  is  based  solely  on  neutron 
counting  rates  above  background  levels.  Gamma  and  X-ray  spectrometers  detect  the 
ionizations  caused  by  those  photons  interacting  within  the  active  detector  volume. 
The  amount  of  ionization  created  by  a  particle  is  dependent  on  its  kinetic  energy 
and  the  energy  required  to  create  charges  within  the  material  through  which  it  is 
passing.  Neutrons  do  not  cause  ionizations  directly  and  as  result,  a  real-time  neutron 
spectrometer  has  not  yet  been  realized  [5].  The  false-positive  and  false-negative  report 
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rates  can  therefore  improve  only  by  more  accurate  and  reliable  neutron  counting-rate 
based  systems. 

An  effort  to  integrate  neutron-sensitive  material  into  solid-state  devices  dates  back 
more  than  twenty  years.  Solid-state  detectors  offer  several  ways  of  improving  neutron 
detector  capabilities  in  areas  where  other  detector  types  cannot.  Scintillation  detec¬ 
tors  have  relatively  poor  energy  resolution  and  slow  timing  characteristics.  Charge 
collection  in  proportional  counters  is  limited  due  to  the  great  differences  in  the  elec¬ 
tron  and  ion  mobilities  in  the  gas.  Ion  mobilities  can  be  orders  of  magnitude  less  than 
electron  mobilities,  and  as  a  result,  the  time  constant  of  detection  circuits  will  only 
allow  for  electrons  to  contribute  to  the  output  signal.  Furthermore,  some  isotopes  im¬ 
portant  in  proportional  counters  designed  to  detect  neutrons,  namely  3He,  are  limited 
in  supply  and  consequentially,  there  is  additional  emphasis  placed  on  finding  alter¬ 
native  neutron  detection  solutions.  Semiconductors  on  the  other  hand,  offer  superior 
energy  resolution,  fast  timing  characteristics,  and  are  relatively  compact.  The  main 
issue  challenging  their  use  as  neutron  detectors  stems  from  the  fact  that  most  popular 
semiconductor  materials  have  low  sensitivity  to  neutrons. 

The  primary  consideration  for  a  detector  is  the  interaction  rate  between  the  de¬ 
tector  medium  and  the  incident  radiation.  This  rate  is  driven  by  the  cross  section 
and  density.  Semiconducting  boron  carbide  has  a  large  thermal  neutron  capture  cross 
section  as  shown  in  Figure  1.  Second  only  to  15'Gd  (considering  only  isotopes  that 
can  form  semiconducting  solids),  the  large  10B  thermal  neutron  capture  cross  section 
improves  the  probability  of  inducing  a  neutron  interaction.  Boron  carbide  is  also  a 
dense  material  (~2.5  [g  cm'3]).  The  combined  effect  of  these  two  qualities  increases 
the  neutron  interaction  rates,  and  reduces  the  neutron  mean  free  paths  in  the  detec¬ 
tor.  The  mean  free  path  is  of  critical  importance  when  developing  compact  detectors 
in  the  sense  that  the  active  volume  of  the  detector  should  be  large  enough  to  stop  a 
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Figure  1.  Total  neutron  cross  section  of  isotopes  common  to  neutron  detection  [6]. 

majority  of  the  particles  incident  upon  it. 

10B  captures  thermal  neutrons  and  emits  'Li  and  4He  particles  according  to  (1). 
It  is  of  no  use  to  have  a  high  probability  of  capturing  neutrons  and  produce  reaction 
products  that  are  still  difficult  to  detect.  This  is  clearly  not  the  case  with  boron;  the 
capture  products  are  highly  energetic  and  can  can  deposit  all  of  their  energy  in  only 
a  few  microns. 

|  94%  4He(1.47  MeV)  +  7Li(0.84  Me V)  +  7(0.48  Me V) 

I  6%  4He(1.78  MeV)  +  7Li(1.02  MeV) 

Furthermore,  boron  has  a  low  Z-number  and  is  not  highly  sensitive  to  gamma  radia¬ 
tion.  Figure  2  shows  that  boron  has  significantly  smaller  cross  sections  in  comparison 
with  other  semiconducting  elements  used  in  radiation  detection.  For  these  reasons, 
semiconducting  boron  carbide  has  been  pursued  for  more  than  fifty  years  [8] .  However, 
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Figure  2.  Photoelectric  and  Compton  Scatter  mass  attenuation  coefficients  for  semi¬ 
conducting  compounds  used  for  neutron  detections  [7]. 


it  is  only  recently  that  material  science  and  technology  capabilities  have  advanced  far 
enough  to  consider  the  employment  of  boron  carbide  in  neutron  detectors. 


1.2  Research  Objectives 

Boron  carbide  devices  have  been  used  in  the  past  to  detect  neutrons.  Most  of 
these  experiments  were  performed  with  boron  carbide  hints  deposited  adjacent  to  the 
active  region  of  a  detector  (conversion-layer  devices  are  discussed  in  Section  2.1.2) 
[9,  10,  11,  12,  13].  Because  the  intrinsic  efficiency  of  these  detectors,  roughly  4%  [14], 
is  low  compared  to  other  detection  methods,  variations  in  geometry  and  fabrication 
techniques  have  been  explored  relentlessly. 

The  diodes  used  in  this  study  were  constructed  with  amorphous  boron  carbide 
films  deposited  on  high  resistivity  n-type  silicon  wafers.  The  objective  in  doing  so 
was  to  create  neutron-sensitive  hetero-j unction  diodes  where  neutron  capture  byprod- 
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ucts  were  created  within  the  active  detection  region.  Due  to  the  complications  and 
challenges  involved  with  amorphous/crystalline  hetero-junction  fabrication,  articles 
claiming  to  have  successfully  detected  neutrons  with  direct-conversion  hetero-junction 
diodes  are  few  in  number.  In  fact,  most  articles  pertaining  to  this  specific  subject  are 
simulations-based  [9]. 

The  objective  of  this  study  was  to  detect  thermal  neutrons  using  amorphous- 
boron-carbide-on-silicon  hetero-junction  diodes.  In  order  to  correctly  assess  the  de¬ 
tectors’  performances,  the  diode  material  and  semiconducting  properties  were  metic¬ 
ulously  studied.  The  ultimate  failure  of  the  diodes  to  detect  neutrons  was  closely 
linked  to  these  properties,  as  well  as  geometric  parameters.  Negative  experimental 
results  usually  lead  to  follow-on  experiments  to  determine  the  root  cause  of  failure. 
In  the  case  of  this  study,  the  causes  for  unsuccessful  neutron  detection  were  pursued 
through  computer  simulation  and  modeling.  Taurus  Davinci  was  used  to  model  the 
device  and  the  electrical  transport  of  injected  charges.  This  simulation  confirmed 
post-characterization  ab  initio  calculations  predicting  that  depletion  widths  were  ad¬ 
equate  to  completely  stop  the  injected  charge,  but  the  weak  electric  fields  and  hence, 
the  slow  drift  velocities,  were  identified  as  the  root  cause  for  poor  charge  collection 
efficiency. 
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II.  Theory 


Neutron  detection  capabilities  have  lagged  behind  detectors  of  other  radiations 
such  as  gamma  or  alpha  particles.  This  is  primarily  due  to  the  lack  of  an  efficient 
solid-state  neutron  detector  [9].  While  still  in  early  stages  of  development,  solid-state 
neutron  detectors  have  been  constructed  from  a  variety  of  elements  capable  of  being 
incorporated  into  semiconducting  solids  and  possessing  high  neutron  cross  sections. 
Most  of  the  successful  diodes  thus  far  have  been  conversion  layer  devices  (Section 
2.1.2).  This  research  investigates  the  use  of  boron  carbide/silicon  hetero- junction 
diodes  to  detect  thermal  neutrons.  The  experimental  approach  of  this  study  stems 
from  fundamental  neutron  detection  principles,  coupled  with  the  physics  of  solid-state 
detectors.  This  chapter  contains  a  brief  review  of  neutron  detection  and  a  discussion 
of  boron  carbide  and  its  properties  that  make  it  suitable  and  appealing  for  neutron 
detection. 

2.1  Neutron  Detection 

Neutrons  carry  no  electronic  charge,  and  therefore  cannot  be  detected  directly 
using  the  measurement  of  electronic  current  or  potential.  Some  neutral  particles,  such 
as  gammas,  interact  with  materials  in  ways  that  cause  ionizations,  and  the  resulting 
electrons  can  then  be  collected  and  measured.  Neutrons  do  not  cause  ionizations 
directly  but  may  result  in  nuclear  reactions.  Some  nuclear  reactions  generate  products 
that  can  ionize  atoms  within  the  detector  volume,  and  this  premise  is  the  basis  for  the 
detectors  used  in  this  research.  Not  all  neutron  interactions  are  useful  in  detection 
however. 

Once  a  neutron  interacts  with  a  target  nucleus,  one  of  two  processes  may  occur. 
The  first  possibility  is  an  elastic  collision.  This  situation  may  change  the  neutron 
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energy  and/or  direction,  but  likely  will  not  prompt  the  production  of  charged  par¬ 
ticles.  In  contrast,  an  atom  can  absorb,  or  capture,  the  incoming  neutron  and  emit 
secondary  radiations  in  the  form  of  heavy  charged  particles.  There  is  a  finite  and 
energy-dependent  probability  of  either  type  of  interaction  to  occur. 

This  probability,  referred  to  as  the  cross  section  per  nucleus,  cr,  is  measured  in  units 
of  barns  (1  b  =  10"24  cm2).  Neutrons  of  a  fixed  energy  have  a  constant  probability 
for  each  nuclear  reaction  mechanism  (such  as  elastic  scattering  or  capture  reactions). 
Because  the  cross  section  may  vary  up  to  several  orders  of  magnitude  depending  on 
the  neutron  energy,  detector  designs  and  methods  differ  greatly  for  neutron  detection 
in  different  energy  regions.  It  is  conventional  to  refer  to  neutrons  with  respect  to 
their  kinetic  energy.  Neutrons  having  kinetic  energies  below  the  cadmium  cutoff  (or 
0.5  eV)  are  designated  as  “slow  neutrons”.  “Fast  neutrons”,  or  neutrons  whose  kinetic 
energies  exceed  0.5  eV,  can  be  detected  with  equipment  and  procedures  different  from 
those  used  for  slow  neutron  detection.  Detection  of  fast  neutrons  is  not  the  focus  of 
this  work,  and  thus  will  only  be  periodically  mentioned. 

Slow  neutrons  scatter  elastically  with  an  atom,  or  experience  one  of  several  neutron 
capture  reactions.  Elastic  scattering  occurs  when  slow  neutrons  collide  with  target 
nuclei  in  the  absorbing  media  but  have  insufficient  kinetic  energy  to  transfer  a  large 
amount  of  energy  to  that  nucleus.  Elastic  collisions  serve  an  important  role  in  neutron 
detection  by  acting  as  a  moderating  force  for  the  neutrons.  This  force  reduces  neutron 
kinetic  energy  into  thermal  equilibrium  with  the  media.  Thermal  energy  is  on  the 
order  of  kT  where  k  is  the  Boltzman  constant  (8.62 xlO"5  eV/K),  T  is  temperature, 
and  correlates  to  0.0258  eV  at  room  temperature  (~300  K). 


2.1.1  Detector  Properties 


Thermal  neutron  detectors  use  materials  with  constituents  that  have  high  cross 
sections  for  thermal  neutron  capture  reactions.  Some  of  the  more  common  isotopes 
incorporated  into  thermal  neutron  detectors  are  shown  in  Figure  1  [6].  Naturally 
occurring  elements  may  actually  be  composed  of  two  or  more  isotopes.  For  neutron 
detector  materials  and  elements,  one  of  the  isotopes  must  have  a  high  thermal  neutron 
capture  cross  section  in  order  to  reliably  lead  to  nuclear  reactions;  the  remaining 
isotopes  in  the  material  may  have  very  low  cross  sections  and  the  detector’s  successful 
operation  is  unaffected  by  their  presence.  Figure  1  shows  the  total  cross  section  for 
specific  isotopes  as  opposed  to  the  cumulative  cross  section  calculated  by  the  ratios  of 
isotopic  abundance  in  naturally  occurring  elements.  The  isotope  pertinent  to  neutron 
detection  must  be  of  high  abundance,  or  an  economical  method  of  enrichment  must 
exist  in  order  for  the  material  to  be  a  viable  means  of  detecting  neutrons. 

In  addition  to  having  large  thermal  neutron  capture  cross  sections,  reactions  em¬ 
ployed  for  slow  neutron  detection  preferably  result  in  heavy  charged  particles  [14]. 
This  is  the  case  with  some  of  the  isotopes  referenced  in  Figure  1;  they  all  have  a 
high  probability  to  capture  neutrons  and  as  indicated  by  (2-6),  these  reactions  all 
have  have  large  Q-values  as  well.  The  Q-value  is  the  amount  of  energy  released  by 
a  reaction.  Reactions  with  large  Q-values  may  deposit  a  greater  amount  of  energy 
within  the  detector  volume. 


^B  +  Jn  3840  b  >  %a  +  1°Li  +  7  +  2.31  MeV  or  2.79  MeV  (2) 

gHe  +  Jn  'j’H  +  )H  +  764  keV  (3) 

^Cd  +  In  “fCd  +  7  +  9.042  MeV  (4) 

jjLi  +  Jn  940  b  >  +  2«  +  4.78  MeV  (5) 

^JG d  +  Jn  Gd  +  7  +  IC_°e  +  7.937  MeV  (6) 
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Neutron  reactive  films  based  on  the  15'Gd(n,7)158Gd  reaction  have  higher  proba¬ 
bilities  than  10B(n,o;) 'Li,  or  6Li(n,a)3H  based  films  of  capturing  neutrons  for  a  fixed 
him  thickness  because  of  the  high  157Gd  thermal  neutron  cross  section.  However, 
the  combined  emission  of  low  energy  gamma  rays,  the  small  branching  ration  of  in¬ 
ternal  conversion  electrons,  and  the  long  range  of  the  internal  conversion  electrons 
make  neutron-induced  events  easily  hidden  by  background  gamma-ray  events.  This 
is  because  only  a  small  number  of  those  internal  conversion  electrons  will  generate 
sufficient  charge  to  exceed  the  intrinsic  noise  level  of  commercial  charge-sensitive 
preamplifiers.  The  particle  energies  emitted  from  the  6Li(n,a)3H  reaction  are  greater 
than  those  emitted  from  the  10B  reaction.  Yet,  the  optimized  him  thickness  for  6LiF 
is  over  ten  times  greater  than  needed  for  10B  with  relatively  unimproved  neutron 
detection  efficiency  due  to  the  lower  cross  section  for  capture  [12]. 

It  is  also  common  to  use  cadmium  foils  as  a  means  of  discriminating  between  those 
signals  induced  by  neutrons  versus  the  signals  caused  by  other  radiation  sources.  As 
will  be  discussed  later,  this  approach  was  used  during  this  research  as  a  means  of 
positively  identifying  neutron  signals. 

One  final  consideration  for  a  detector  is  its  sensitivity  to  gamma  radiation.  Gam¬ 
mas  are  ubiquitous  to  most  neutron-rich  environments  and  therefore  must  be  ac¬ 
counted  for.  As  shown  in  (2),  (4),  and  (6),  gammas  can  be  a  product  of  the  reaction 
ultimately  responsible  for  the  generating  a  signal.  Gamma  rays  can  ionize  and  ex¬ 
cite  atoms  within  the  detector;  both  cases  result  in  unwanted  detector  output  which 
can  ultimately  complicate  the  task  of  detecting  neutrons.  Neutron  detectors  should 
therefore  be  constructed  of  material  with  low  gamma  sensitivity,  but  must  also  have 
a  geometry  that  will  also  prevent  large-scale  gamma  interference.  For  example,  most 
neutron  proportional  counter  tubes  are  smaller  in  diameter  than  the  mean  free  path 
of  most  gamma  rays.  Gammas  can  then  pass  through  the  detector  and  deposit  only 
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a  portion,  if  any,  of  their  energy. 


2.1.2  Detector  Technologies 

Radiation  detectors  typically  fall  into  one  of  three  classes:  proportional  counters, 
scintillators,  and  solid-state  detectors.  Although  the  term  solid-state  is  somewhat 
ambiguous  since  scintillation  detectors  could  technically  be  considered  solids,  the 
accepted  practice  is  to  categorize  only  those  detectors  whose  signal  is  based  on  the 
collection  of  electron-hole  pairs  as  solid-state  detectors.  More  simply  stated,  solid- 
state  devices  are  almost  always  semiconducting  detectors,  and  for  this  reason,  both 
terms  will  be  used  interchangeably  throughout  this  document. 

The  purpose  of  this  research  is  to  improve  existing  methods  for  detecting  neu¬ 
trons.  An  understanding  of  earlier  methods  used  in  neutron  detection  work  will  help 
gauge  the  success  of  current  research  efforts,  and  to  what  degree.  Therefore,  this 
section  contains  of  brief  overview  of  non-solid  state  neutron  detectors  in  addition  to 
a  more  in-depth  explanation  of  this  work’s  central  theme,  semiconductor  properties 
and  detectors. 

2. 1.2.1  Non-solid-state  detectors 

Proportional  counters  used  in  neutron  detection  are  commonly  cylindrical  tubes 
filled  with  3He  or  BF3  gas.  Figure  3  depicts  the  basic  proportional  counter  config¬ 
uration.  An  electric  held  is  established  in  the  detector  by  applying  high  voltage  to 
the  anode.  When  exposed  to  thermal  neutrons,  charged  particles  created  by  (2)  or 
(3)  ionize  gas  molecules  in  the  detector,  resulting  in  electron  drift  towards  the  an¬ 
ode.  Strong  electric  fields  can  generate  a  cascade  of  ionizations  when  electrons  gain 
sufficient  kinetic  energy  to  cause  ionizations  when  they  collide  with  gas  molecules. 
The  secondary  electrons  are  accelerated  by  the  electric  held.  This  gas  multiplication 
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Anode 


Figure  3.  Basic  elements  of  a  proportional  counter.  The  cathode  ensures  there  is  a 
vacuum-tight  enclosure  to  contain  the  fill  gas.  Positive  voltage  is  applied  to  the  anode 
wire  which  creates  an  electric  field  perpendicular  to  the  longitudinal  axis  (z-axis)  [14]. 

process  is  known  as  the  Townsend  avalanche. 

Scintillators  are  commonly  used  in  fast  neutron  detection,  but  also  can  be  used  in 
slow  neutron  detection  with  inorganic  lithium-containing  crystals.  The  crystal  lattice 
in  the  scintillator  determines  the  energy  states  of  the  material.  Although  electrons 
are  typically  found  in  the  valence  band,  neutron  absorptions  can  excite  electrons 
into  the  conduction  band.  The  products  of  the  6Li(n,T)4He  reaction  (5)  deposit 
4.78  MeV  to  the  crystal  lattice.  Photons  are  emitted  by  the  process  of  fluorescence 
when  the  electrons  relax  from  the  conduction  band  to  the  valence  band.  Photo¬ 
multiplier  tubes  are  commonly  incorporated  into  scintillator  detectors  in  order  to 
collect  photons  and  convert  them  via  photo-sensitive  layers  into  photoelectrons  which 
can  then  be  multiplied  in  order  to  yield  a  measurable  current  pulse.  Although  beyond 
the  scope  of  this  research,  scintillators  are  commonly  used  in  fast  neutron  detection 
as  well.  This  form  of  detection  relies  on  the  recoil  protons  produced  when  high-energy 
neutrons  collide  elastically  with  atoms  in  the  scintillator  lattice. 
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2. 1.2. 2  Solid-state  detectors 


Solid-state,  or  semiconductor,  detectors  rely  on  the  formation  of  electron-hole 
pairs,  analogous  to  electron-ion  pairs  in  gas-filled  proportional  counters,  to  generate 
a  signal.  Ionizing  radiation  interacts  with  solid  matter  in  the  detector  and  deposits 
energy,  elevating  valence  band  electrons  to  the  the  conduction  band.  Physically,  this 
simply  represents  an  electron  gaining  sufficient  energy  to  leave  the  bonding  site  and 
drift,  under  an  electric  held,  through  the  crystal  lattice  [14].  This  process  creates 
a  vacancy,  or  hole,  in  the  valence  band.  The  amount  of  energy  required  to  create 
an  electron-hole  pair,  or  the  ionization  energy,  is  dependent  on  the  band  gap  of  the 
semiconductor.  Temperature  is  another  factor  in  the  semiconductor  ionization  energy, 
but  to  a  much  smaller  degree.  With  an  applied  electric  held,  both  the  electrons  in 
the  conduction  band  (holes  in  the  valence  band)  will  respond  to  electrostatic  forces 
and  result  in  a  net  migration  of  charge. 

The  motion  of  both  electrons  and  holes  will  be  a  combination  of  both  diffusion 
and  drift  velocity  [14].  The  motion  of  electrons  will  be  anti-parallel  to  the  electric 
held,  but  the  drift  velocity  of  holes  is  conceptually  more  difficult  to  grasp.  The 
hole  represents  the  absence  of  a  negative  charge,  thereby  creating  a  net  positive 
charge.  Both  electron  and  hole  drift  velocities  are  proportional  to  the  electric  held  at 
low  to  medium  held  strengths.  With  sufficiently  large  electric  helds,  the  saturation 
velocity  is  reached.  Section  3.3.1  will  examine  this  topic  more  closely,  but  charge 
collection  time  corresponds  to  the  time  required  for  the  charges  to  transport  from  the 
point  of  generation  to  the  depiction  region  boundary.  Most  detectors  are  operated  at 
saturation  so  that  the  charged  particles  will  be  accelerated  from  the  depletion  region 
with  the  greatest  possible  velocity,  thereby  minimizing  the  charge  collection  times 
(ideally  on  the  order  of  a  ns). 

Electron  or  hole  movement  constitutes  a  current  that  persists  until  the  charge 
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carrier  reaches  the  boundary  of  the  detector’s  active  volume.  Assuming  energy  is 
deposited  at  a  single  point,  the  resulting  carrier  currents  will  still  be  different  be¬ 
cause  the  charge  collection  times  vary  with  the  carrier  mobilities.  In  some  modern 
semiconductor  detectors,  hole  mobilities  are  close  enough  to  electron  mobilities  such 
that  both  carrier  currents  can  be  integrated  on  a  measuring  circuit  with  a  long  time 
constant.  This  is  not  feasible  in  proportional  counters  where  the  ion  mobility  is  orders 
of  magnitude  greater  than  the  electron  mobility  [14]. 

Solid-state  detectors  emphasized  in  this  work  are  based  on  the  properties  near  a 
junction  of  n-  and  p-  type  semiconducting  material.  The  effects  of  carrier  diffusion 
across  this  p-n  junction  gives  rise  to  an  electric  held  and  a  region  depleted  of  majority 
charge  carriers,  the  depletion  region.  The  premise  of  diode  radiation  detectors  is  that 
the  electron-hole  pairs  created  within  the  depletion  region  are  swiftly  accelerated  out 
of  the  region  by  the  electric  held.  The  motion  of  the  charges  produces  an  output 
current  pulse.  This  transient  current  is  measurable  only  if  the  steady  state  leakage 
current  is  sufficiently  low.  Ideally,  the  leakage  currents  should  be  on  the  order  of  a 
nA,  or  10"9  A  [14].  The  pulse  rise  time  is  defined  by  migration  time  from  the  point 
of  formation  to  the  boundary  of  the  depiction  width. 

In  the  absence  of  an  applied  voltage,  the  natural  built-in  voltage  across  the  de¬ 
pletion  width  causes  the  carriers  to  drift,  but  the  performance  of  the  detector  will 
be  inherently  poor  under  these  circumstances.  This  is  due  to  the  higher  probability 
for  recombination  and  charge  trapping  when  the  carriers  move  slowly  through  the 
depletion  region,  an  event  the  depletes  the  charge  collection  potential. 

The  fraction  of  charges  escaping  capture  increases  as  the  electric  held  becomes 
stronger.  Reverse  biasing  the  diode  increases  the  the  electric  held  and  the  thickness 
of  the  depletion  region  [15].  This  depletion  region  represents  the  active  volume  of  the 
detector.  As  will  be  discussed  in  Section  3.3.2,  devices  used  in  this  study  were  only 
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Figure  4.  Electric  field  comparison  between  fully-  and  partially-depleted  detectors  [16] 


partially  depleted.  Most  detectors  are  operated  fully-depleted  for  several  reasons. 
Figure  4  effectively  illustrates  the  differences  in  the  electric  fields  for  detectors  that 
are  partially  depleted  (shown  on  the  left)  versus  those  that  are  fully  depleted  (shown 
on  the  right).  In  the  figure,  w  signifies  the  edge  of  the  depletion  boundary  and 
corresponds  to  the  distance  from  the  junction  where  the  electric  held  vanishes,  which 
has  deleterious  effects  on  the  drift  velocity,  pulse  rise  time,  and  the  increases  the 
probability  of  charge  trapping.  The  un-depleted  dead  layer  between  the  depletion 
region  and  the  electrode  is  also  source  of  noise  and  can  deteriorate  the  detector 
resolution. 

In  fully-depleted  detectors,  the  electric  held  exists  throughout  the  diode  as  a  whole 
because  depletion  region  spans  the  whole  width  between  the  two  electrodes.  As  the 
applied  bias  increases,  the  electric  held  increases,  and  the  difference  between  Emax 
and  Emin  is  minimized.  Eventually,  if  the  bias  continues  to  increase,  the  electric  held 
will  become  uniform  throughout  the  diode  [16]. 

Solid-state  thermal  neutron  detectors.  Solid-state  thermal  neutron  de¬ 
tectors  incorporate  material  with  high  cross  sections  for  neutron  capture  interactions. 
Conversion  layers  and  solid-state  direct  conversion  devices  (Figure  5)  both  employ 
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Figure  5.  Comparison  of  (a)  conversion  layer  and  (b)  direct  conversion  diodes  [5]. 

neutron  induced  reactions.  The  most  common  reactions  used  to  produce  detectable 
ionization  are  the  10B(n,a)'Li,  6Li(n,a)3H,  and  the  113Cd(n,7)114Cd  reactions  [12]. 

Planar  conversion  layer  devices  use  neutron-sensitive  thin- films  deposited  adjacent 
to  a  diode  under  reverse  bias.  Reaction  products  emitted  close  enough  to  and  in 
the  direction  of  the  diode  junction  can  create  electron-hole  pair  production  within 
the  depletion  region  of  the  reverse-biased  diode.  The  electric  field  present  in  the 
detector  minimizes  the  recombination  of  the  electron-hole  pairs.  As  the  thin-hlm 
thickness  increases,  so  does  the  probability  of  capturing  a  neutron  and  producing 
charged  particles.  These  reaction  products  travel  finite  distances  in  material  prior 
to  depositing  all  of  their  energy.  Thin-hlms  that  are  too  thick  may  absorb  all  of  the 
charged  particle  energy  before  the  particle  reaches  the  detector  volume.  An  additional 
factor  concerning  conversion  layer  devices  is  the  inherent  cap  in  intrinsic  efficiency. 
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The  incoming  thermal  neutron  energy  (0.0259  eV)  is  quite  low  compared  to  the  Q- 
value  of  the  reaction.  Both  resultant  products  are  emitted  in  opposite  directions  in 
order  to  conserve  momentum.  Thus,  there  is  a  high  likelihood  that  only  one  reaction 
product  will  be  emitted  in  the  direction  of  the  detector  volume. 

Direct  conversion,  or  solid- form  [12],  devices  overcome  the  disadvantages  charac¬ 
terizing  conversion  layer  devices  by  integrating  the  neutron-sensitive  thin  film  within 
the  diode  junction.  In  semiconductor  detectors,  the  measurable  signal  originates  from 
electron-hole  pair  formation  in  the  depletion  region.  The  obvious  advantage  to  solid 
form  devices  is  that  the  neutron  interaction  occurs  within  the  active  detector  region 
and  all  reaction  products  are  capable  of  producing  detectable  ionizations.  Factors 
such  as  the  film  thickness,  the  location  and  direction  of  the  emitted  particles,  bias 
voltages,  recombination  rates,  and  the  material  properties  of  the  film  itself  all  deter¬ 
mine  the  probability  of  producing  a  measurable  signal. 

2.2  Properties  of  Boron  Carbide 

Few  elements  possess  a  neutron  capture  cross  section  suitable  for  neutron  de¬ 
tectors.  Even  fewer  elements  meet  this  criteria  and  are  also  capable  of  forming 
semiconducting  solids.  Boron  accomplishes  both  criteria.  10B  is  20%  abundant  in 
naturally  occurring  boron  and  has  a  large  capture  cross  section  for  thermal  neutrons. 
Boron  also  forms  a  boron-rich  semiconducting  solid,  boron  carbide  [17],  which  exists 
in  many  forms.  Most  of  the  literature  addressing  the  semiconducting  properties  of 
boron  carbide  focus  on  the  crystalline  form  with  B:C  ratios  ranging  from  4.3:1  to  11:1 
[18,  19,  20,  21,  22,  23].  The  boron  carbide  incorporated  into  the  devices  used  in  this 
study  were  amorphous,  which  is  the  least  well  referenced  form  of  boron  carbide.  This 
gap  in  knowledge  precludes  direct  comparisons  between  devices  cited  in  the  literature 
and  the  devices  used  in  this  study  without  taking  into  account  how  the  amorphous 
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Figure  6.  Typical  rhombohedral  lattice  structure  in  semiconducting  boron  carbide  [24]. 


structure  affects  the  physical  properties. 

While  amorphous  structures  lack  a  recurrent  long-range  order,  the  semiconductor 
properties  of  the  material  do  not  appear  to  depend  upon  crystallite  size  [8].  Sim¬ 
ilarities  between  crystalline  and  amorphous  boron  carbide,  as  will  be  reported  in 
Section  2.2.1,  will  be  basis  of  the  following  discussion.  The  significance  of  the  simi¬ 
larities  leads  to  a  logical  conclusion:  The  principles  associated  with  crystalline  boron 
carbide,  such  as  crystal  structure  and  electrical  transport  phenomena,  are  also  closely 
related  to  the  principles  of  the  amorphous  counterpart. 

2.2.1  Structure 

Crystalline  BC  structures  vary  but  are  often  icosahedral-crystal  type  structures 
with  rhombohedran  lattice  structures  as  shown  in  Figure  6.  The  crystal  structures 
of  boron  carbides,  including  amorphous  structures,  share  the  icosahedron  building 
block  that  is  also  found  in  rhombohedral  boron  [25].  The  icosahedral  is  the  pri¬ 
mary  conduction  mechanism  for  the  material.  The  presence  of  two  carbon  atoms, 
serving  as  electron  donors,  slightly  distorts  the  molecular  polyhedral  cage  struc¬ 
ture  [26].  Figure  7  shows  the  molecular  structure  of  ortho-  and  meta-carborane 
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(a) 


Figure  7.  Molecular  structure  of  (a)  ortho-  and  (b)  meta-carborane  C2B10H12.  Carbon 
atoms  are  shaded  [26]. 


atoms.  When  either  molecule  is  used  as  a  chemical  vapor  deposition  precursor,  the 
chiral  differences  between  the  isomers  lead  to  different  majority-carrier  semiconduc¬ 
tors.  Meta-  carborane  (closo-l,7-dicarbadodecaborane)  forms  an  n-type  semiconduc¬ 
tor  while  ortho-carborane  (closo-l,2-dicarbadodecaborane)  forms  a  slightly  p-type 
semiconductor  [27].  Ortho-carborane  deposited  by  Plasma-Enhanced  Chemical  Va¬ 
por  Deposition  (PECVD),  forms  an  amorphous  semiconductor  denoted  as  a-E>5C:Hx. 
A  schematic  picture  of  this  structure  is  shown  in  Figure  8. 

The  semiconducting  properties  of  boron  carbide  make  it  uniquely  suitable  for  de¬ 
tection  diodes.  Semiconducting  films  can  be  incorporated  in  a  p-n  junction,  thereby 
removing  the  requirement  for  a  neutron  sensitive  conversion  layer  adjacent  to  the 
diode  [5,  9].  Furthermore,  the  consequence  of  forming  different  majority-carrier 
semiconductors  can  be  exploited  to  create  an  all-boron  carbide  heteroisomeric  diode. 
Meta-carborane/ortho-carborane  used  as  deposition  precursors  in  sequence  have  pro¬ 
duced  rectifying  diodes  [29]. 
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Figure  8.  Schematic  depiction  of  a— B5C:HX  [28]. 


2.2.2  Boron  Carbide  Electrical  Transport 

Boron  carbide  forms  micro-crystalline  structures,  and  no  unit  cell  describes  the 
structure  as  a  whole  [18].  The  electrical  transport  properties  of  boron  carbide  are 
not  well  understood,  but  one  theory,  from  the  works  of  Werheit  et.al..,  is  the  current 
dominant  theory.  Preceding  Werheits  theory  is  an  explanation  based  on  a  model 
proposed  by  Emin  et.al.  [30].  From  his  works,  Emin  argues  that  bipolaron  hopping 
is  responsible  for  boron  carbide  electrical  activity.  Recent  experimental  evidence 
published  by  Werheit  et.al.  contradict  the  theory  and  provide  findings  incompatible 
with  the  presence  of  bipolarons  in  boron  carbide  [19,  21,  22,  23]. 

Werheits  model  contends  the  semiconducting  nature  is  derivative  from  the  high 
density  of  anti-site  and  vacancy  defects  [20].  The  theory  follows  from  theoretical 
band  structure  calculations  and  the  proposed  band  diagram  is  shown  in  Figure  9  [23]. 
Point  defects  in  semiconductors  generate  split-off  valence  states  in  the  band  gap. 
These  unoccupied  valence  states,  according  to  calculated  electron  deficiencies  and  ex¬ 
perimentally  determined  point  defect  densities,  are  exactly  compensated  in  real  solids. 
This  implies  both  that  the  material  is  a  semiconducting  solid  and  the  valence  bands 


20 


r-  0 


2 


<L> 


LU 

>N  1 

05 

<D 

c 

LU 


0 


CoofluconPand 


1  56 01 
1  5695 

1  2 

0  92 
077 

047 


electron  traps 
•  excitons 


0  27 


-  0.5 
f  E(eV) 


Range 

of 

deep  levels 


max.  total  density: 
*1  /  elan,  cell 


0.18 

0.065 


-3-10'W 


Valence  bam 


t 


- (acceptor  level) 


Figure  9.  Energy  band  scheme  based  on  optical  absorption,  luminescence,  XRS,  and 
transport  properties.  Left  ordinate,  energies  related  to  the  valence  band  edge;  right 
ordinate,  energies  relative  to  the  conduction  band  edge.  Arrows  indicate  the  direction 
of  the  measured  optical  transition  [23]. 


are  completely  filled.  The  electron  deficiency  evokes  the  generation  of  compensating 
intrinsic  defects,  likely  for  energetic  reasons.  Such  high  concentrations  of  gap  states 
attached  to  the  valence  band  affect  the  electronic  charge  transport;  essentially,  they 
are  responsible  for  the  p-type  character  and  the  very  low  electrical  conductivity  [18] . 

Dowben  et.al.  [31]  have  attempted  to  determine  the  band  gap  of  boron  carbide 
via  optical  absorption  measurements  at  varying  ratios  of  boron  to  carbon  (Figure  10). 
The  results  indicate  the  band  gap  is  closely  correlated  with  the  boron  to  carbon  ratio 
and  can  range  from  0.77  to  1.80  eV  and  is  consistent  with  the  thermal  activation 
barrier  of  1.25  eV  for  conductivity  [31].  This  conclusion  differs  from  the  2.09  eV 
band  gap  of  Figure  9.  Any  direct  comparison  between  the  reports  however  is  possibly 
flawed.  First,  it  is  possible  the  0.9  eV  band  gap  indicated  by  the  optical  absorption 
measurements  is  a  non-direct  transition  between  the  gap  states  and  conduction  band. 
No  attempt  is  made  in  [31]  to  clarify  whether  the  reported  values  reference  non-direct 
transitions  or  indirect  transitions  from  the  valence  to  the  conduction  band.  Several 
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Figure  10.  Band  gaps  determined  (for  a  BsC:Hx)  at  room  temperature  are  plotted 
against  composition  determined  by  Auger  electron  spectroscopy.  The  measurement 
are  largely  determined  by  optical  absorption  •  but  data  for  photoluminescence  A  is 
also  shown  [31]. 

non-direct  transitions  determined  through  optical  absorption  in  polycrystalline  boron 
carbide,  including  one  0.92  eV  transition,  are  mentioned  in  [23]. 

High  electron  deficiency  induced  gap  states  results  in  high  carrier  concentration  on 
the  order  of  1018  -  1019  as  shown  in  Figure  11(a)  [20].  Drude-type  transport  coupled 
with  a  hole  hopping  mechanism  within  the  partially  filled  gap  states  is  responsible 
for  transporting  charge  carriers  through  the  boron  carbide  [21]. 

The  large  number  of  gap  states  increases  the  probability  of  trapping.  As  shown  in 
Figure  11(b),  Werheit  reports  boron  carbide  ambient  mobilities  converging  to  values 
on  the  order  of  1  [cm2/V-s]  [20]. 

However,  charge  collection  measurement  times  typical  of  p-type  a-B5C:Hx  thin 
films  on  n-type  silicon  diodes  exposed  to  a  neutron  held  result  in  calculated  mobilities 
on  the  order  of  10"2  -  10"4  [cm2/V-s]  [17].  The  reduced  mobility  in  thin  films  is  likely 
due  to  greater  defect  site  concentrations  expected  of  amorphous  solids  versus  poly- 
or  mono-crystalline  solids.  Table  1  summarizes  some  of  the  electronic  properties 
discussed  herein. 
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Figure  11.  (a)  Comparison  of  carrier  densities  and  (b)  Carrier  mobility  in  boron  carbide 
as  reported  by  different  experimental  methods  and  authors  [20]. 


2.3  Previous  Work 


Boron  carbide  has  been  shown  to  form  excellent  lietero-,  and  hetero-isomeric  junc¬ 
tions  that  have  stable  leakage  currents  close  to  ideal  levels  [32],  Some  of  the  neutron 
detection  studies  performed  in  the  past  give  great  insight  in  how  to  conduct  the 
current  research,  as  well  as  interpret  and  analyze  data. 


Table  1.  Electronic  properties  of  boron  carbide 


Parameter 

Symbol 

Value 

Werheit 

Dowben 

Mobility 

Total  Carrier  Density 
Dopant  Concentration 
Band  gap  energy 

fip 

l^n 

Nt 

Na 

Eg 

1  cuP/V-s  [18] 

1018  -  1019  cm"3 

2.09  eV  [23] 

10"2  -  10"4  cnr/V-s  [17] 
0.001  cm2/V-s  [17] 

[18] 

4.5 xlO12  cm"3  [17] 

0.77  -  1.88  eV  [31] 
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Conversion-layer  devices.  Conversion- layer  devices  were  among  the  first 
approaches  to  implement  semiconductor  devices  for  neutron  detection.  The  neutron 
capture  products  are  only  produced  in  the  the  neutron-sensitive  him;  for  the  pur¬ 
poses  of  this  discussion  the  him  is  boron  carbide.  As  the  him  thickness  increases, 
the  probability  of  capturing  incident  neutrons  increases  and  approaches  unity  with 
sufficiently  thick  conversion-layer  thickness.  In  other  words,  the  capture  efficiency  is 
theoretically  100%. 

While  the  neutron  capture  efficiency  improves  with  the  thicker  films,  a  point  is 
reached  where  the  overall  detector  efficiency  suffers  as  result  of  excessive  him  thick¬ 
nesses.  The  'Li  and  4 He  particles  deposit  energy  immediately  from  the  point  of 
generation.  All  of  the  energy  deposited  by  the  charged  particles  in  the  boron  carbide 
region  will  be  lost  due  to  recombination.  If  the  particles  deposit  all  of  their  energy 
by  this  manner,  they  are  said  to  have  “self  absorbed.”  Thus,  it  is  possible  to  cap¬ 
ture  thermal  neutrons  with  great  efficiency  but  still  have  none  of  the  resulting  energy 
deposited  in  the  detector. 

In  contrast,  particles  born  within  a  infinitesimally  thin  him  will  likely  deposit  only 
a  small  fraction  of  their  energy  before  reaching  the  detector.  These  competing  factors 
pose  a  difficult  dilemma  in  determining  the  best  detector  design  using  conversion-layer 
films.  McGregor  et.al.  have  explored  this  problem  in  great  detail  with  numerical 
models  and  boron  carbide  conversion-layer  device  simulations  and  as  result,  have 
reported  optimized  parameters  regarding  the  specihc  geometry  of  devices. 

Figure  12(a)  and  12(b)  summarize  their  findings.  Figure  12(a)  shows  that  the 
least  energetic  particle,  'Li  at  840  keV,  deposits  all  of  its  energy  in  less  than  2  /mi. 
In  comparison,  the  predominant  a  particle  of  the  reaction  will  deposit  870  keV 
(1470  keV  -  600  keV)  over  the  same  distance  [11].  Figure  12(b)  represents  the  opti¬ 
mized  parameters  that  factor  in  to  the  overall  detector  efficiency.  Due  to  the  differ- 
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Figure  12.  (a)7Li  and  4He  range  (in  boron  carbide)  and,  (b)  film  thickness  optimization 
determined  from  a  computational  analysis  [11]. 


ences  in  track  length  shown  in  Figure  12(a),  the  optimum  hint  thickness  varies  with 
respect  to  each  of  the  four  particles  [11].  With  that  said  however,  it  is  safe  to  conclude 
that  the  intrinsic  efficiency  suffers  with  hints  greater  than  2.5  /irn. 

The  pulse-height  spectra  taken  front  conversion- layer  diodes  typically  exhibit  two 
or  more  spectral  features.  A  perfect  example  of  this  is  presented  in  Figure  13  [33]. 
Both  reaction  products  are  emitted  back  to  back,  and  only  one,  if  either,  forms 
electron-hole  pairs  within  the  detector.  The  two  dominant  peaks  correspond  to  'Li 
and  4He  at  840  keV  and  1.47  MeV.  The  1.02  MeV  7Li  is  likely  masked  within  the 
continuum  of  the  4 He  peak  at  1.47  MeV. 


Direct  conversion  devices.  The  feature  most  notable  regarding  conversion- 
type  devices  is  the  absence  of  a  full-energy  peak,  or  a  peak  corresponding  to  the  total 
Q- value  of  the  reaction.  Most  of  the  work  done  thus  far  dealing  with  direct-conversion 
diodes  has  been  computer  model-based.  In  fact,  many  authors  claiming  to  have  suc¬ 
cessfully  detected  neutrons  via  direct-conversion  devices  have  incorrectly  analyzed 
pulse-height  spectra  data  by  failing  to  note  the  absence  of  the  full  energy  peak  [12]. 
Figure  14  [9]  shows  that  the  the  diode  geometries  play  a  significant  role  in  the  detec- 
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Figure  13.  Pulse  height  spectra  collected  using  conversion-layer  diode  [33]. 

tor  response  to  neutrons.  Immediately,  one  notices  two  significant  differences.  First, 
both  spectra  exhibit  a  continna  around  the  Li  and  He  peaks,  but  the  conversion  layer 
signatures  have  tailing  towards  lower  energies  and  direct-conversion  devices  have  tails 
that  extend  towards  higher  energies.  In  conversion  layer  devices,  the  particles  deposit 
energy  before  reaching  the  detector  volume,  and  the  particle  generation  site,  direc¬ 
tion,  and  path  length  all  contribute  to  an  energy  continuum  below  the  full  energy 
value  of  the  particle.  A  similar  effect  is  observed  in  direct  conversion  devices.  This 
is  explained  by  the  following  argument.  Whenever  charged  particles  are  created,  one 
particle  will  likely  deposit  only  some  of  its  energy  while  the  other  makes  a  full  en¬ 
ergy  deposition.  The  energy  from  both  particles  is  collected  simultaneously,  but  path 
length  of  the  particle  will  create  a  continuum  as  well.  In  this  case  however,  the  energy 
is  integrated  with  the  energy  of  the  oppositely  directed  particle. 

Secondly,  the  direct-conversion  device  shows  two  full  escape  peaks  for  both  re¬ 
action  possibilities.  These  peaks  cannot  occur  in  conversion  devices.  As  the  him 
thicknesses  grow  larger,  the  spectra  tend  to  flatten  at  lower  energies,  and  the  two 
full  energy  peaks  grow  in  count  number.  In  other  words,  the  thickness  of  the  him  is 


26 


0.01 


Figure  14.  Computer  simulation  prediction  of  both  a  conversion-layer  and  direct  con¬ 
version  diode  spectra  [9]. 


sufficiently  thick  to  stop  both  particles  within  the  active  detector  volume. 
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III.  Experimental  Approach 


3.1  Device  Fabrication 

This  study  was  performed  using  three  different  devices.  Two  of  devices  were 
hetero- junction  diodes.  One  diode  was  fabricated  at  the  University  of  Missouri, 
Kansas  City  (UMKC).  The  second  diode  was  fabricated  in  2010  at  the  University 
of  Nebraska,  Lincoln  (UNL).  It  is  worth  noting  that  LINL  diode  was  used  in  conjunc¬ 
tion  with  thesis  work  in  2010.  A  resistive  type  device  was  the  second  class  of  device 
tested  during  this  research;  it  was  also  constructed  at  UNL.  The  underlying  physics 
and  materials  science  involving  the  fabrication  processes  for  these  devices  is  beyond 
the  scope  of  this  thesis,  but  a  brief  description  of  the  fabrication  process  is  included. 
This  chapter  comments  briefly  on  the  fabrication  procedures  for  the  devices  used  in 
this  research  (a  more  detailed  procedure  can  be  obtained  from  [34,  35]),  and  followed 
by  the  device  precharacterization  and  experimental  procedures. 

3.1.1  Hetero-junction  Device 

UMKC  Diode.  The  diodes  were  prepared  in  a  custom-built  capacitively 
coupled  PECVD  reactor  by  the  reaction  of  ortho-carborane  vapor  in  a  200  mTorr 
argon  background.  The  use  of  PECVD  provides  the  means  for  fabricating  boron 
carbide  thin  films  at  lower  temperatures  and  pressures  (<1000  °C,  <50  Torr)  [31].  The 
boron  carbide  is  assumed  to  have  a  much  greater  carrier  concentration  than  the  high 
resistivity  silicon.  The  dopant  concentration,  Np,  for  10  [kO-crri]  n-type  silicon  is  on 
the  order  of  1011  cm"3  [15].  The  function  of  high  resistivity  (low  dopant  concentration) 
silicon  is  two-fold.  The  overall  space  charge  neutrality  must  be  maintained  (the  total 
space  charge  on  either  side  of  the  junction  is  equal),  and  is  given  by  (7).  xp  and  xn 
correspond  to  the  depletion  thickness  on  the  p-  and  n-side  of  the  junction,  respectively. 


The  greater  carrier  concentration  on  the  boron  carbide  side,  NAl  will  deplete  a  greater 
volume  on  the  high  purity  silicon  side. 


Naxp  =  NDxn  (7) 

Dopants  can  also  act  as  traps  within  the  active  detector  volume.  Thus,  high  pu¬ 
rity  silicon  results  in  greater  depletion  widths  for  relatively  thin  boron  carbide  films. 
Charge  carriers  produced  as  a  result  of  energetic  charged  particles  will  be  less  likely 
to  become  captured  by  traps  within  the  active  detector  volume. 

High  resistivity  (>10  [kD-cm])  n-type  silicon  (111)  wafers  with  a  thickness  of 
515±20  /im  were  used  for  the  substrate.  The  wafers  were  diced  into  square  pieces 
of  15x15  mm  for  film  deposition  and  device  fabrication  purposes.  Before  the  vapor 
deposition  process  could  be  performed,  the  substrates  were  cleaned  using  a  Piranha 
solution  to  remove  the  organic  impurities  from  the  surface  and  were  chemically  etched 
in  hydrofluoric  acid  to  remove  surface  oxides.  The  cleaned  substrates  were  rinsed  with 
acidified  deoxygenated  water  and  dried  using  IPA  vapor.  This  five-step  sequence 
ensured  organic  impurities,  particulates  and  surface  oxides  were  minimized. 

Once  cleaned  and  dried,  the  substrates  were  inserted  into  a  carousel  cradle  as¬ 
sembly,  capable  of  holding  nine  substrates.  The  assembly  rotated  during  deposition 
for  film  uniformity.  Care  was  taken  to  ensure  all  nine  substrates  were  loaded  into 
the  plate  assembly  and  transferred  to  the  primary  PECVD  reactor  chamber  within 
15  minutes  to  prevent  new  oxides  from  developing  on  the  clean  surfaces.  The  cradle 
assembly  formed  the  positive  electrode  in  the  PECVD  system  and  was  powered  by 
RF  generator  operating  at  13.56  MHz.  Using  a  custom  solid-source  bubbler,  solid 
ortho-carborane  precursor  was  sublimed  at  72  -  75  °C  and  delivered  to  the  reactor 
with  an  argon  carrier  gas.  The  gaseous  argon/carborane  mixture  was  released  into 
the  reactor  chamber  through  a  shower  head,  which  also  formed  the  negative  electrode 
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to  the  capacitively  coupled  PECVD  system  [34], 

Once  the  boron  carbide  films  were  grown,  elipsometric  measurements  were  per¬ 
formed  to  determine  the  him  thickness.  The  him  had  a  thickness  of  1.5±0.2  pm.  A 
secondary  sputtering  chamber  was  used  to  deposit  aluminum  contacts  on  the  diode. 
The  contacts  on  the  boron  carbide  side  were  deposited  using  a  0.2  cm2  shadow  mask. 
The  boron  carbide-Al  layer  was  protected  using  photoresist  and  the  reverse  side  of 
the  substrate  was  cleaned  using  hydrofluoric  acid  solution,  and  was  rinsed  and  dried 
using  the  procedure  above.  This  was  done  to  ensure  the  surface  oxides  had  been 
removed  prior  to  contact  deposition.  The  photoresist  was  removed  by  dipping  the 
device  in  acetone  solution  after  the  contacts  were  deposited  on  the  reverse  side. 

The  device,  now  having  aluminum  contacts  on  both  the  boron  carbide  and  silicon 
side,  could  be  characterized  and  inspected  for  measured  to  determine  if  it  functioned 
as  a  diode.  It  was  first  mounted  on  a  fiberglass  board  with  two  separate  copper-clad 
inlays  (shown  in  Figure  15).  The  diode  was  affixed  to  the  largest  of  the  two  copper- 
clad  inlays.  Conductive  silver  paint,  lightly  dabbed  on  the  silicon  (n-side)  contact  of 
the  diode,  served  as  both  a  means  of  connectivity  between  the  n-side  of  the  diode  and 
the  copper-clad  contact,  as  well  as  a  bonding  agent  to  secure  the  diode  in  place.  A 
thin  copper  wire  connected  the  boron  carbide  (p-side)  contact  to  the  smaller  copper 
inlay. 


UNL  Diode.  Specifics  regarding  the  geometry  of  UNL  diode  are  estimated. 
The  diode  differs  from  UMKC  diode  in  geometry,  most  notably  in  the  boron  carbide 
thickness  and  the  contact  composition  and  size.  The  diode  was  created  via  PECVD; 
a  summarized  discussion  of  this  process  can  be  found  earlier  in  this  section.  Hemiseal 
was  used  to  coat  the  device  which  prevented  him  thickness  measurements.  All  other 
differences  between  UMKC  and  UNL  diodes  can  be  found  in  the  Table  2.  The  LINL 
diode  was  secured  to  the  fiberglass  board  in  accordance  with  the  procedure  described 
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Figure  15.  Close-up  photo  of  the  diode.  Aluminum  contacts  were  sputtered  onto  both 
sides,  and  was  mounted  on  fiberglass  board.  Silver  paint  secured  the  the  n-side  contact 
to  the  largest  copper  inlay.  A  thin  copper  wire  connected  the  boron  carbide  (p-side) 
to  the  smaller  copper  inlay.  The  BNC  center-pin  wire  is  shown  soldered  to  the  p- 
side  contact.  The  wire  was  detached  and  soldered  to  the  n-side  inlay  for  the  neutron 
irradiation  phase  of  this  study. 


for  UMKC  diode. 


Table  2.  Diode  Material/ Geometric  Parameters 


UMKC  Diode 

UNL  Diode 

Silicon  Resistivity 

>10  [kO-crn] 

8  [kO-crn]  [36] 

Silicon  Thickness 

515  [/irn] 

525  H 

BC  Resistivity 

108  -  1010  [0-ern]  [17] 

108  -  1010  [Q-crn]  [17] 

BC  Thickness 

1.5  [/mi] 

2.0  [/mi] 

Contact  Area 

2.0  xlO'1  [cm2] 

3.32  xlO'2  [cm2] 

Contact  Metal 

Aluminum 

Silver  chromate 

3.1.2  UNL  Resistive- type  Device 

The  UNL  resistor-like  device  incorporates  the  neutron  capturing  potentials  of 
157Gd,  6Li,  and  10B.  Films  of  Gd203  and  Li2B407  were  deposited  on  sapphire  glass 
(A1203)  via  pulsed  laser  deposition  (PLD)  [35].  The  gadolinium  oxide  film  was  limited 
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to  60  nm  in  order  to  avoid  excessive  gamma  interference.  Li2B407  was  deposited  on 
the  opposite  side  of  the  device.  The  PLD  procedure  was  not  a  vital  part  of  this 
study,  the  method  of  deposition  is  detailed  in  [35].  As  with  the  diodes,  a  copper-clad 
fiberglass  board  was  used  to  secure  the  device.  Two  1.3  x  10"2  cm2  gold  contacts  on  the 
gadolinium  side  of  the  him  (facing  away  from  the  fiber-glass  board)  were  successfully 
bound  to  the  copper-clad  firmly  using  conductive  silver  paint. 

3.2  Pre-Irradiation  Device  Characterization 

The  Keithley  4200  semiconductor  characterization  system  performed  all  prechar¬ 
acterization  measurements  during  the  course  of  this  research  via  the  the  Signatone 
CM-220  probing  station.  All  measurements  were  performed  at  approximately  300  K. 

3.2.1  Diode  I— V  and  C— V  characterization 

I— V  The  p-n  junction  serves  as  the  rectifying  element  in  the  diode,  allowing 
current  to  how  freely  in  one  direction  and  inhibiting  it  in  the  other.  Current-voltage 
(I-V)  measurements  were  performed  on  both  diodes  to  verify  rectifying  behavior  and 
to  characterize  the  diode  leakage  current.  Due  to  the  device  area  differences,  a  more 
useful  comparison  was  the  current  density  (J).  In  Figure  16,  both  diodes  exhibit 
rectifying  behavior.  The  y-axis  (anode  current  density)  is  set  on  a  logarithmic  scale. 
The  error  bars  are  smaller  than  the  pixel  diameter  of  the  plotted  data.  It  is  apparent 
that  differences  in  time,  fabrication  location,  and  diode  geometry  affected  the  I-V 
characteristics,  manifested  by  the  differences  in  the  current  densities  for  all  applied 
anode  voltages.  Figure  16  shows  that  the  the  forward  and  reverse  bias  currents 
densities  differ  by  three  orders  of  magnitude  for  the  UMKC  diode.  The  UNL  diode 
rectification  was  considerably  smaller  with  only  two  orders  of  magnitude  difference 
between  the  forward  and  reverse  bias  currents. 
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Figure  16.  J— V  curve  for  UMKC(blue)  and  UNL(green)  diodes  (y-axis  is  a  log  scale). 
The  predicted  Js  for  the  UMKC  diode  is  plotted  in  black. 

Both  diodes  exhibited  stable  currents  between  1.0  -  5.0  V  reverse  bias.  The  UNL 
diode  produced  a  current  density  of  1.49(±0.03)xl0~6  and  1.24(±0.03) xlO'4  [A/cm2] 
at  5.0  V  reverse  and  forward  bias  respectively.  The  corresponding  values  for  the 
UMKC  diode  are  1.07(±0.06) xlO"4  and  1.83(±0.2) xlO"1  [A/cm2].  These  values  were 
compared  to  the  the  theoretical  reverse  saturation  current  density  in  order  to  quantify 
the  leakage  current. 

Total  diode  current  is  given  by  the  Shockley  equation,  also  known  as  the  ideal 
diode  law  (8).  The  equation,  as  written,  was  modified  to  account  for  different  intrinsic 
carrier  densities,  rq,  of  the  different  semiconductor  materials  involved  in  the  hetero¬ 
structures  [37].  D  and  L  are  the  diffusion  constants  and  lengths  respectively  for  the 
minority  carriers.  Na  and  Njj  are  the  accepter  and  donor  concentrations  in  boron 
carbide  and  silicon  respectively.  Under  forward  bias,  the  current  density  increases  at 
a  constant  rate,  and  under  reverse  bias  however,  the  current  saturates  at  the  reverse 
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saturation  current. 


J 


Dnn\bc)q  |  DPn\si)q 
LhNa  LpNd 


Js 


(8) 


Leakage  current,  caused  by  thermionic  emission  and  carrier  tunneling,  has  adverse 
effect  on  the  signal-to-noise  (S/N)  ratio  when  using  diodes  as  radiation  detectors. 
In  reference  to  Section  2.1.2,  leakage  currents  on  the  order  of  nA  are  ideal.  The 
saturation  current  density,  Js,  was  computed  from  values  listed  in  Table  4  and  is 
plotted  in  black  in  Figure  16.  Appendix  B.1.1  also  contains  this  calculation.  Some 
parameters  determined  experimentally  did  not  agree  with  values  reported  in  literature 
sources.  Other  parameters  had  large  uncertainties  due  to  the  wide  range  of  possible 
values.  In  these  cases  the  mid-range  value  was  used  for  the  approximation  (e.g. 
effective  hole  mass  in  boron  carbide).  Precise  leakage  current  calculations  could  not 
be  obtained  for  the  UMKC  diode  due  to  these  uncertainties.  The  leakage  current  was 
estimated  to  be  ~10  /j,A  at  -5  V,  and  calculated  from  (8)  with  the  assumption  that 
excess  conductivity  above  the  reverse  saturation  current  was  due  to  leakage. 

When  the  experimental  data  are  plotted  against  the  Shockley  equation  (8),  there 
is  poor  agreement  between  the  measured  values  and  the  theoretical  predictions  of 
the  forward  current  as  shown  in  Figure  17(a).  The  ideal  diode  equation  assumes  an 
ideality  factor  (r/)  of  1.0.  This  is  consistent  with  forward  I-V  characteristics  when 
diffusion  current  dominants;  q  increases  as  recombination  current  becomes  dominant. 
A  more  accurate  expression  for  diode  current  density  is 


J 


Js 


(9) 


A  least  squares  fit,  shown  in  purple  in  Figure  17(b),  of  the  linear  forward  current 
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corresponds  to  77=1.0  (the  forward  current  is  linear  only  with  a  semi-log  scale).  If  we 
extrapolate  the  line  towards  the  y-axis,  the  intercept  occurs  at  the  reverse  saturation 
current  density  ( J(0V)  =  Js).  The  relatively  weak  correlation  between  the  theory  and 
experimental  data  is  caused  by  a  larger  ideality  factor  for  the  UMKC  diode.  Using 
the  the  same  y-intercept,  a  line  tangent  to  the  measured  data  was  projected  (shown 
in  green  in  Figure  17(b)).  The  slope  of  the  line  corresponded  to  ideality  factor  of 
3.5.  The  forward  current  sub-linearity  prevented  an  effort  to  calculate  the  ideality 
factor  empirically,  or  via  the  method  described  using  the  ideal  diode  forward  current 
linear  region.  Shockley-Read-Hall  recombination  theory  predicts  an  ideality  factor 
equal  to  2.0  or  less.  The  deviation  from  this  prediction  is  likely  attributed  to  the 
high  defect  states  of  the  amorphous  boron  carbide.  The  high  defect  density  of  boron 
carbide,  on  the  order  of  1021  cm'3  [20],  makes  multi-level  recombination  possible.  In 
the  case  of  high  defect  densities,  ideality  factors  can  increase  above  Shockley-Read- 
Hall  recombination  predictions  of  1.0  <  77  <  2.0  [38]. 

Given  the  qualitative  ideality  factor  analysis,  it  is  not  surprising  that  both  diodes 
exhibited  sub-linear  forward  current  behavior.  This  is  consistent  with  a  high  resis¬ 
tivity  material.  Boron  carbide  diodes  characterized  in  [32]  were  compared  to  other 
hetero-junction  diodes  (ZnSe/ZnCdSe)  and  the  sub-linearity  of  the  forward  current 
was  interpreted  as  an  effective  resistor  in  series  (both  the  amorphous  boron  carbide 
and  n-type  silicon  have  large  impedences).  The  boron  carbide  films  are  highly  resistive 
(on  the  order  of  109  Q-crri  at  room  temperature),  and  the  forward  current-limiting 
behavior  shown  in  Figure  16  can  be  attributed  the  diodes’  highly  resistive  amorphous 
region  [31,  32], 

Due  to  the  shortage  of  devices  exhibiting  proper  rectifying  I-V  characteristics, 
break-down  voltage  analysis  was  omitted.  A  diode  reverse  biased  to  breakdown  can 
show  increases  in  leakage  current.  In  diode  detectors,  depleting  the  diode  to  the 
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Figure  17.  (a)  Current  density  of  UMKC  diode  vs  Shockley  ideal  hetero-structure 

model,  (b)  UMKC  diode  vs  Shockley  ideal  model  with  linear  region  fits  accounting  for 
ideality  factors. 


greatest  possible  extent  is  beneficial  because  it  increases  the  active  detector  volume 
and  improves  charge  collection  efficiency.  Determining  the  optimal  detector  bias  may 
result  in  damaging  diodes  and  given  the  limited  number  of  samples  available  for 
testing  during  this  study,  breakdown  analysis  could  not  be  performed. 
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C— V  Under  reverse  bias,  the  diode  junction  acts  like  a  parallel-plate  ca¬ 
pacitor.  In  order  to  calculate  the  depletion  width  of  the  diode,  capacitance-voltage 
measurements  (C-V),  measurements  were  performed.  Capacitance  can  be  related  to 
the  depletion  width  of  a  diode  via  (10)  where  Cj  is  the  junction  capacitance,  A  is  the 
diode  area,  es  is  the  semiconductor  permittivity,  and  W  is  the  depletion  width. 


Cj  =  — 

J  w 


(10) 


The  C-V  measurements  plotted  in  Figures  18(a)  and  18(b)  for  each  diode  indicate 
the  neither  device’s  C-V  characteristics  are  ideal.  The  left  ordinate  corresponds  to  the 
measured  capacitance  (converted  to  units  of  F  cm"2);  the  right  ordinate  corresponds 
to  the  inverse  squared  data  (1/C2).  The  UMKC  diode  capacitance  exceeds  the  UNL 
diode  capacitance  by  1  -  2  orders  of  magnitude  depending  on  the  applied  reverse  volt¬ 
age.  Depletion  inversely  affects  the  junction  capacitance  (10),  and  the  large  UMKC 
diode  capacitances  indicate  the  diode  is  less  depicted  than  than  the  UNL  diode.  Fur¬ 
thermore,  achieving  the  smallest  detector  capacitance  reduces  electronic  noise  and 
can  ultimately  improve  the  detector  resolution  [14].  Thus,  detectors  are  operated 
with  the  largest  possible  depiction  region  in  order  to  maximize  the  radiation-sensitive 
volume,  and  to  improve  the  S/N  ratio  by  minimizing  the  capacitance. 

Junction  capacitance  can  be  rewritten  in  terms  of  the  semiconductor  permittiv¬ 
ities,  dopant  concentrations,  and  the  applied  voltage  (11).  The  1  / \]Vu  —  V  depen¬ 
dency  corresponds  to  asymptotic  increase  in  Cj  as  the  applied  reverse  bias  approaches 
0  V.  It  can  also  be  inferred  from  (11)  that  junction  capacitance  does  not  exist  when 
the  forward  bias  exceeds  V/,; .  From  this  perspective,  it  is  apparent  the  UNL  diode  has 
non  diode-like  C-V  characteristics.  In  Figure  18(b),  the  capacitance  increases  as  the 
reverse  bias  decreases  (as  expected).  However,  the  peaking,  followed  by  a  decrease  of 
the  capacitance  as  the  applied  voltage  approaches  zero  is  indicative  of  some  phenom- 
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(a) 


(b) 

Figure  18.  Capacitance  measurements  for  each  diode.  The  left/right  ordinate  refers  to 
the  C— V / C'2  measurement  for  (a)  UMKC  diode,  and  (b)  UNL  diode. 


ena  other  than  junction  depletion.  Although  the  direct  cause  of  this  behavior  was  not 
investigated,  it  is  possibly  linked  to  radiation  damage  effects  incurred  from  previous 
research. 


Cj  = 


qesiebcNDNA 


(11) 


2  (csiNd  +  ^bcNA)  (Vbi  —  V) 

Equation  11  can  be  manipulated  algebraically  so  the  function  is  linear  with  respect 
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Figure  19.  l/Cj2  of  UMKC  diode(blue)  with  linear  fit  (green). 

to  the  independent  variable,  V.  The  inverse-squared  of  (12)  for  an  abrupt  hetero¬ 
junction  diode  results  in  a  straight  line.  Figures  18(a)  and  18(b)  differ  in  this  aspect. 
Only  the  UMKC  diode  1  /Cj2  plotted  data  is  consistent  with  an  abrupt  junction;  the 
corresponding  UNL  diode  data  produces  a  curved  line  which  contradicts  the  abrupt 
junction  assumption.  The  fabrication  process  for  the  UNL  diode,  discussed  within 
[36],  was  the  basis  for  this  conclusion. 


a 


2  {(-siNo  +  (-bc^A)  (Hi  —  V) 
qesiebcNDNA 


(12) 


Figure  19  presents  the  UMKC  diode  1  /Cj2  data  with  a  linear  fit.  The  expression 
for  this  fit  is  included  in  the  upper  right  corner  of  the  figure.  The  correlation  between 
the  the  two  lines  further  supports  that  the  junction  is  abrupt.  A  linear  fit  to  the 
UNL  diode  1  / Cj 2  data  was  not  possible  without  excessive  statistical  uncertainty.  Hi 
can  be  calculated  from  x-axis  intercept  of  the  fitted  line.  Using  (12),  when  V—Vbll 
1/Cj2=  0.  Appendix  B.1.1  provides  a  detailed  calculation  for  determining  the  UMKC 
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diode  built-in  voltage.  Unfortunately,  the  lack  of  linearity  exhibited  by  the  UNL  diode 
1  /Cj2  data  precluded  a  similar  fit  for  the  built-in  voltage.  In  order  to  estimate  the 
individual  depletion  widths  into  the  BC  (where  the  doping  density  is  not  well  known) 
and  the  silicon,  the  BC  region  was  initially  assumed  to  be  fully  depicted  (a  reasonable 
assumption  based  on  [17]).  The  dopant  concentration  was  much  higher  in  the  BC 
region  than  in  the  silicon  [20],  and  the  slope  of  the  1  /Cj2  versus  V  line  was  purely  a 
result  of  depletion  expanding  in  the  silicon. 


(13) 


Cj  for  an  abrupt  hetero-junction  is  treated  as  the  equivalent  capacitance  of  a 
two-capacitor-in-series  connection.  The  equivalent  capacitance  (13)  was  used  to  ap¬ 
proximate  the  hetero-junction  depletion  width  [39].  Using  the  measured  capacitances 
(converted  to  F  cm"2),  and  the  parameters  listed  in  Table  3,  the  depletion  width  into 
the  silicon,  xn,  was  calculated. 


Table  3.  Variables  used  in  depletion  width  calculation 


Parameter 

Symbol 

Value 

Permittivity  of  Free  Space 

e0 

8.85  xlO"14  F  cm"1 

Silicon  permittivity 

11.9  e0 

BC  permittivity 

^bc 

8.0  e0  [40] 

Junction  capacitance 

Cj 

(Measured)  F  cm'2 

BC  Thickness 

Xp 

1.5  (UMKC)  /irn 

2.0  (UNL)  //rn 

Applying  reverse  bias  to  a  diode  widens  the  depletion  width  further  into  the  sil¬ 
icon  towards  the  n-side  electrode,  and  increases  the  internal  electric  held  within  the 
depletion  region.  The  detector  efficiency  benefits  from  both  consequences  of  reverse 
bias.  The  result  of  the  depletion  width  calculations  is  shown  in  Figures  20(a)  and 
20(b).  Tables  4  and  5  summarize  the  calculations  derived  from  the  capacitance  mea¬ 
surements  for  both  diodes.  The  active  detector  volume  predictions,  the  UMKC  diode 
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(a)  UMKC  diode  depletion 


(b)  UNL  diode  depletion 


Figure  20.  Depletion  width  comparison  under  reverse  bias. (a)  UMKC.  (b)  UNL 


reverse  saturation  current  density  approximation,  and  the  diode  modeling  during  the 
final  stages  of  this  study  were  based  on  these  values  listed  in  the  two  tables.  N&  and 
Nd  were  only  calculated  for  the  UMKC  diode.  The  experimentally  determined  Np 
exceeded  dopant  concentrations  expected  for  >10  [kf2-cm]  n-type  silicon  as  reported 
by  [15]  by  2  orders  of  magnitude.  The  boron  carbide  effective  dopant  concentrations 
reported  by  [17]  coincided  with  the  UMKC  diode  Np  calculation  within  a  factor 
of  1.6. 


41 


Table  4.  UMKC  Diode  Characterization  Summary 

Parameter  Value 
Vbi  [V]  0.23852 

Nd  [cm'3]  3.13X1013 

Na  [cm-3]  7.0  xlO12 

Wov  [/mi]  10.9402 

W_5v  H  44.5261 

Table  5.  UNL  Diode  Characterization  Summary 

Parameter  Value 
Vbi  [V]  UNK 
Na,d  [cm'3]  UNK 
Wov  t/H  128.804 
W_5V  [/an]  136.073 

The  depiction  region  volume  is  given  by  the  product  of  the  diode  area  and  the 
depletion  width.  The  active  detector  volumes  are  8.9 xlO"4  cm3  and  4.5xl0"4  cm3  for 
the  UMKC  and  UNL  diode  respectively  at  5  V  reverse  bias. 

While  the  advantages  of  greater  depletion  width,  namely  a  stronger  electric  field 
and  greater  detector  volume,  are  clear,  the  device  must  be  depleted  sufficiently  so  that 
ionizing  radiation  can  deposit  most  of  its  kinetic  energy  within  the  active  detector 
volume.  The  a  particles  are  the  most  energetic  reaction  products,  and  they  will  also 
have  longest  ionization  track  in  both  the  boron  carbide  and  silicon.  The  stopping 
power  in  silicon  limits  the  track  length  to  6.0  microns.  The  depletion  widths  in  both 
diodes  exceed  this  length,  and  therefore,  the  detector  is  depleted  adequately  and  will 
allow  the  charged  particle  products  to  deposit  most,  if  not  all,  of  their  kinetic  energies. 

3.2.2  Resistor-like  device  characterization 

The  UNL  resistive-detector  was  ultimately  the  only  device  to  produce  spectra 
with  identifiable  features.  The  I-V  characteristics  of  the  device,  shown  in  Figure  21, 
are  inconsistent  to  a  degree  with  Ohm’s  Law,  V  =  IR.  In  the  case  of  a  resistor, 
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the  expected  I-V  dependency  is  linear.  The  data  plotted  in  Figure  21  was  collected 
by  sweeping  the  voltage  from  -20  -  20  V.  A  slightly  different  curve  resulted  from  a 
reverse  sweep  (higher  to  lower  voltage). 

The  observed  hysteresis  vanished  at  -5  V,  and  as  voltage  was  increased,  the  current 
gradually  became  linear  with  voltage.  The  inverse  effect  was  observed  for  the  high  to 
low  voltage  sweep  (not  shown).  In  this  case,  hysteresis  vanished  at  +5  V.  Plotting 
both  data  sets  together  showed  both  linear  regions  were  in  fact  parallel,  with  a  slope 
coefficient  that  corresponded  to  R"1.  This  demonstrated  that  the  device  was  purely 
resistive,  and  for  the  purpose  of  clarity,  only  one  of  the  sweeps  was  included  in  the 
figure  to  represent  the  I-V  characteristics  without  needlessly  convoluting  the  figure. 

Based  on  the  slope  of  the  linear  region,  the  resistance  of  the  device  is  8.93  xlO10  Q. 
An  issue  left  unresolved  is  the  slight  conductivity  at  0  V.  This  was  partially  due 
to  charge  trapping  and  surface  transport  [41].  The  pre-characterization,  performed 
using  a  probe  station,  was  not  performed  in  complete  darkness.  Aluminum  oxide  is  a 
transparent  material,  and  when  exposed  to  light,  could  have  also  contributed  to  this 


Figure  21.  I— V  measurements  resistor  like  device  with  linear  fit  of  the  I— V  (linear 
region  existed  for  Vappiie(j  greater  than  -15  V). 
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departure  from  the  expected  resistor  characteristics. 


3.3  Irradiation  Experimental  Procedure 

The  goal  of  the  irradiation  experiment  was  to  determine  if  a  hetero-junction  diode 
constructed  of  a-B5C:Hx  and  n-type  silicon  could  be  used  to  detect  thermal  neutrons. 
Theory  and  modeling  supports  that  the  devices  used  for  this  study  were  properly 
configured  to  measure  the  products  of  the  10B  neutron  capture  reaction  [5,  36].  These 
devices  were  fundamentally  different  from  conversion  type  detectors  and  both  of  the 
reaction  products  ('Li  and  a  particles)  should  have  ionized  carriers  in  the  active 
detector  region  in  spite  of  being  emitted  180°  from  each  other. 

Many  modifications  were  made  until  reaching  the  electrical  connections  and  com¬ 
ponent  configuration  described  below.  Only  the  details  pertinent  to  this  final  con¬ 
figuration  are  discussed.  The  over-arching  concept  behind  the  experiment  was  that 
any  neutron- induced  spectral  feature  observed  in  the  unshielded  pulse-height  mea¬ 
surement  would  be  suppressed  with  a  cadmium  foil.  Cadmium  is  a  highly  efficient 
absorber  of  thermal  neutrons  (Figure  1),  and  with  the  device  shielded  with  the  foil, 
only  energetic  neutrons  reach  the  detector.  The  boron  carbide  cross  section  was  low 
for  fast  neutrons,  and  thus  the  probability  of  capturing  the  neutron  was  considerably 
less. 

3.3.1  Electrical  connections 

The  devices  were  packaged  in  an  aluminum  case  as  shown  in  Figure  22.  The  cases 
were  constructed  from  3-in  diameter  aluminum  rods  by  creating  a  small  central  cavity 
inside  to  secure  and  protect  the  device  during  experiments.  A  small  hole  was  drilled 
and  tapped  in  the  side  in  order  to  accept  a  threaded  BNC  connector.  The  purpose  of 
this  design  was  three-fold:  it  prevented  damage  to  the  device,  muted  electrical  noise, 
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Figure  22.  Device  mounted  in  aluminum  case. 


and  formed  a  large  electrical  ground.  The  BNC  center  coaxial  pin  was  insulated  from 
the  aluminum  case,  but  the  ground  (the  outer  metal  wall  of  the  connector)  formed 
a  direct  (tension)  contact  with  the  case  once  it  was  securely  tightened  into  the  side- 
wall  of  the  aluminum  case.  Four  holes  were  drilled  and  tapped  in  the  bottom  of  the 
case’s  internal  cavity  to  match  the  spacings  of  the  holes  drilled  through  the  fiberglass 
boards.  The  two  holes  drilled  through  corners  with  the  copper-clad  inlays  served  an 
additional  purpose.  The  contact  could  be  grounded  to  the  case  by  tightening  a  steel 
bolt  through  the  hole,  which  established  a  pressure  contact  between  the  copper  inlay 
and  the  aluminum  case  (Figure  22).  Nylon  offsets  screwed  into  the  holes  formed  a 
mounting  platform  for  the  fiberglass  boards.  The  nylon  also  functioned  as  an  insulator 
and  ensured  no  part  of  the  board  came  in  contact  with  the  case. 

Figure  23  depicts  the  electrical  connections  for  the  experiment.  A  small  gauge 
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Figure  23.  Electronic  connectivity  setup. 


copper  wire  was  used  to  connect  the  coaxial  center-pin  to  the  n-side  copper  contact. 
All  spectral  data  was  collected  from  the  silicon  (n-side)  with  the  BC  (p-side)  grounded. 
Neutron  capture  occurred  only  in  the  p-side  of  the  device;  however,  most  of  the  charge 
was  captured  in  the  silicon  [17].  Electron-hole  pairs  were  created  from  the  ionizing 
energy  loss  of  the  highly  energetic  'Li  and  4He  particles  and  the  predominant  signal 
collected  at  the  electrode  was  due  to  ionizations  in  the  depletion  region.  Accordingly, 
electrons  (holes)  were  accelerated  towards  the  positive  (negative)  electrode. 

Electron  mobility,  /in  in  silicon  is  1350  [cm2/V-s],  whereas  hole  mobility  in  boron 
carbide  is  much  less  (approximately  0.00075  [cm2/V-s])  [17].  The  vast  differences  in 
carrier  mobility  also  results  in  significant  differences  in  the  charge  collection  times. 
As  the  holes  transport  through  the  boron  carbide,  their  low  mobility  enhances  charge 
trapping  probability.  Charge  trapping  and  de-trapping  broaden  the  distribution  of 
(and  protract)  the  times  required  for  holes  to  reach  the  negative  electrode. 

Charge  carrier  motion  within  the  depletion  region  causes  the  dominant  detector 
output  signal.  The  largest  fraction  of  the  depletion  region  extends  into  the  silicon. 
The  high  purity  of  the  silicon  also  enhances  the  probability  of  charge  carriers  avoiding 
capture  by  traps  and  recombination  centers.  These  two  factors  improve  the  charge 
collection  capabilities  of  the  positive  electrode,  and  justify  collecting  all  spectra  from 
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the  n-side. 


It  is  necessary  for  the  charge  collection  time  (16)  to  be  smaller  than  signal  pro¬ 
cessing  time  constant  (Section  3.3.2)  for  the  signal  pnlse  to  be  a  faithful  measure  of 
the  energy  deposited  [14].  The  charge  collection  times,  tcouection,  were  approximated 
by  (14  -  16). 


£(x) 

QNd  (ur  X 

=  - (IT  —  x) 

(14) 

Vn{x) 

=  Hn£  (x) 

(15) 

[0  1 

t collection 

=  /  — ——ax 

Jx o  Vn[X) 

log 


W 


qHnNd  W  -  x0 


(16) 


The  electric  held  £  during  partially-depleted  detector  operation  is  not  constant.  The 
maximum  electric  held  occurs  at  the  p-n  junction  and  then  decreases  linearly  on  both 
sides  of  the  junction  through  the  depletion  region.  The  held  vanishes  at  the  depletion 
width  boundaries.  (16)  assumes  mobility  is  constant,  and  states  that  the  collection 
time  for  charge  created  at  Xo  is  dependent  on  the  drift  velocity,  vn,  which  also  varies 
with  x  (distance  from  the  depletion  boundary).  The  values  for  electron  mobility  in 
Si,  nn,  and  hole  mobility  in  boron  carbide,  /ip  are  reported  earlier  in  this  section  and 
in  Table  1  respectively. 

The  charge  collection  time  can  also  have  significant  effects  on  the  transient  current 
measured  in  the  detector.  The  charge  collection  time  represents  the  increment  of  time 
over  which  the  deposited  charge  is  integrated.  The  transient  current,  calculated  from 
(17),  must  exceed  the  steady  state  leakage  current  in  order  to  have  a  measurable  pulse 
and  good  S/N  ratios.  The  energy  deposited,  Edep,  is  determined  from  the  reaction 
Q-value,  and  Epri  is  the  energy  required  to  produce  electron-hole  pairs  (3.6  eV  in 
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silicon).  ttrans  is  calculated  from  (16),  and  tpiasma  is  time  needed  to  disperse  the  cloud 
of  charged  pairs  such  that  they  do  not  self-shield  (2-5  ns)  [14,  36].  Calculations 
shown  in  Appendix  B.2  indicate  that  hole  collection  time  is  six  orders  of  magnitude 
greater  than  the  electron  collection  time. 
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3.3.2  Signal  processing  component  configuration 

Before  spectra  were  collected,  the  electronic  connections  and  configurations  of  the 
detection  electronics  shown  in  Figure  24  were  confirmed  with  a  pulse  generator.  The 
pulser  was  connected  to  the  oscilloscope  and  adjusted  to  deliver  a  700  mV  pulse.  Each 
component  was  tested  using  the  pulse  generator  signal.  The  final  step  in  validating 
the  equipment  setup  was  to  observe  the  pulser  spectrum.  The  peak  width  (Figure  25) 
was  only  two  channels,  and  centered  just  below  channel  70.  Channel  70  corresponded 
to  an  input  voltage  into  the  multi-channel  buffer  (with  the  conversion  gain  set  to 
1024  channels)  of  700  mV.  This  agreed  with  the  700  mV  oscilloscope  pulse,  and  the 
fine  width  of  the  peak  indicated  that  connections  were  correct,  and  that  no  observ¬ 
able  electronic  component-induced  noise  entered  the  detection  circuitry.  The  specific 
settings  and  other  pertinent  information  in  reference  to  the  electronic  modules  were 
listed  in  Appendix  A. 

Several  preamplifiers  were  considered  for  this  experiment  and  it  is  worthwhile 
discussing  the  factors  leading  to  the  final  selection.  The  Ortec  142IH  preamplifier 
is  a  charge  sensitive  universal  purpose  instrument  that  can  be  used  for  many  classes 
of  detectors  including  semiconductor  detectors  [42],  The  instrument  also  accepts 
input  from  a  bias  supply  to  apply  voltage  to  the  detector.  By  this  configuration,  the 
diode  bias  is  sent  to  the  detector  and  the  resulting  response  signal  is  delivered  to  the 


48 


Figure  24.  Experimental  procedure  equipment  connection  diagram 
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Figure  25.  Output  spectrum  for  the  pulser  producing  700  mV  pulse  on  oscilloscope. 


preamplifier  input  over  the  same  cable. 

The  preamplifier  time  constant,  r,  was  the  determining  factor  in  this  experi¬ 
ment.  The  time  constant  is  the  product  of  the  input  impedance  and  dynamic  capaci¬ 
tance,  which  for  the  142IH  are  93  and  10,000  pF  respectively  [42],  The  case  where 
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(r  >>  tc)  ensures  that  pulse  generated  is  proportional  to  the  charge  deposited  in  the 
detector  [14].  Appendix  B.2  contains  the  values  of  tc  and  r  with  the  appropriate 
calculations.  Calculations  in  Appendix  B.2  indicate  collection  times  on  the  BC  side 
were  greater  than  the  time  constant. 

The  background  spectra  were  collected  for  each  detector  prior  to  placing  it  near 
the  source.  Spectra  were  then  measured  for  300/600/900  s  while  being  irradiated 
by  the  neutron  source.  A  low  level  discriminator  was  used  to  ensure  dead  time  did 
not  exceed  10%.  Once  complete,  a  cadmium  foil  was  wrapped  around  the  aluminum 
case  and  an  additional  300  s  spectra  was  collected.  The  pulse-height  spectra  from 
UMKC  diode  showed  spectral  features  with  these  spectra  duration  times  but  due  to 
the  low  energy  of  the  pulses,  the  interference  from  background  was  significant.  An 
additional  set  of  spectra  was  collected  for  UMKC  diode  only,  but  the  duration  times 
of  the  spectra  were  increased  to  600/1800/3600  s. 

3.3.3  Neutron  source 

The  neutron  irradiation  experiment  was  performed  at  The  Ohio  State  University 
Nuclear  Reactor  Lab.  A  5.0  Ci  Plutonium-Beryllium  (PuBe)  source  was  used  as  a 
source  of  neutrons.  The  PuBe  source  contained  a  mixture  of  Pu  and  Be.  The  source 
was  sealed  in  an  15  gallon  steel  drum  filled  with  borated-parrafin.  A  central  cavity 
of  the  “neutron  howitzer”  was  hollow  and  lined  with  a  7  5/s”  diameter  aluminum 
pipe.  The  PuBe  source  was  stored  inside  that  aluminum  pipe,  and  no  moderating 
material  was  placed  above  the  source  in  the  original  configuration  [43].  Two  1-in  discs 
of  polyethylene  were  placed  between  the  source  and  the  detector  while  spectra  were 
collected. 

Plutonium  decays  and  produces  alpha  particles  and  uranium  (18  -  19).  Beryllium 
absorbs  those  a’s  and  produces  carbon  and  neutrons  (22).  Additional  neutrons  born 
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from  (21)  increase  the  neutron  emission  rate  for  several  years  after  the  production  of 
the  source  [44], 
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Use  of  the  Bateman  equations  is  required  to  calculate  the  neutron  emission  rate  for 
(18)  thru  (22).  The  neutron  emission  rate  calculation,  modified  from  [45],  is  included 
in  Appendix  B.3.  Figure  26(a)  shows  that  the  PuBe  neutron  spectrum  is  predomi¬ 
nately  fast  neutrons,  with  the  largest  peak  centered  at  4  MeV  [46].  The  curve  maxi- 
mas  in  Figure  26(b)  correlate  to  the  thickness  of  polyethylene  required  to  therm  alize 
a  beam  of  mono-energetic  neutrons  [47].  The  dashed  green  corresponds  to  1  MeV 
fast  neutrons,  which  become  thermalized  after  penetrating  3.0  cm  of  polyethylene. 
By  interpolating,  4.0  cm  of  polyethylene  (1.08  cm  less  than  the  thickness  of  the  two 
polyethylene  discs  used  during  the  experiment)  is  an  ample  thickness  for  moderating 
most  of  the  PuBe  source  neutron  population. 

The  geometry  of  the  experiment  is  shown  in  Figure  27.  Neutrons  emitted  were 
captured  by  the  borated-paraffin  unless  they  were  emitted  vertically  upward  through 
the  void  in  the  parraffin  directly  above  the  source.  The  detector  was  therefore  placed 
8  inches  above  the  source.  Activity  was  calculated  at  1.19x10'  [n  s'1].  With  the 
known  isotropic  neutron  emission  rate,  the  thermal  flux  at  the  detector  location  as 
well  as  the  rate  reaction  density  could  be  calculated  (included  in  Appendix  B.3). 
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(a)  (b) 

Figure  26.  (a)PuBe  neutron  spectrum  [46].  (b)  Thermalization  as  a  function  of 

polyethylene  thickness  [47]. 
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Figure  27.  Geometry  of  the  irradiation  experiment  using  a  PuBe  source.  The  source 
and  detector  were  separated  by  eight  inches  (including  two  inches  of  polyethylene). 
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The  thermal  flux  was  estimated  using  (23),  where  A  is  the  neutron  emission  rate 
and  B  is  the  branching  ratio  for  the  reaction  of  interest.  From  (1),  the  branching  ratio 
was  94%.  0/47 r  is  the  probability  that  the  source  neutron  hits  the  detector,  and  12 
is  the  the  solid  angle  subtended  by  the  detector  [48].  The  thermal  macroscopic  cross 
section,  S,  is  given  by  (24),  where  ath  is  the  microscopic  thermal  capture  cross  section 
and  N  is  the  atomic  density.  Given  a  thermal  flux  of  464  [n/cm2-s],  the  reaction  rate 
density  was  then  computed  as  1.14x10'  [cm-3  s'1]. 
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IV.  Results  and  Analysis 


4.1  Irradiation  Results 

None  of  the  devices  used  in  this  study  produced  spectral  features  or  oscilloscope 
traces  that  could  be  linked  to  neutron  capture  events.  The  UNL  resistive  device  spec¬ 
tra  had  some  features,  but  the  increased  pulse  amplitudes  and  frequencies  observed  in 
the  oscilloscope  while  the  device  was  shielded  by  a  cadmium  foil  conclusively  indicated 
that  neutron  capture  products  were  not  producing  the  observed  signal. 

The  UMKC  diode  did  not  produce  spectra  with  identifiable  features.  Pulses  ob¬ 
served  on  the  oscilloscope  indicated  the  possibility  of  an  interaction,  but  when  com¬ 
pared  to  the  pulse  height  spectra,  it  was  apparent  the  pulses  were  due  mostly  to 
noise  and  gamma/X-ray  interference.  Recall  that  the  UMKC  diode  exhibited  diode¬ 
like  C-V  characteristics  but  pre-characterization  measurements  also  indicated  high 
leakage  currents  (~10  /i A  at  5.0  V  reverse  bias)  and  large  capacitances.  These  two 
undesirable  properties  of  the  diode  likely  contributed  to  the  failure  to  detect  thermal 
neutrons. 

The  UNL  diode  did  not  yield  a  single  useable  spectrum.  The  300  s  measurement 
totaled  more  counts  in  several  of  the  low  energy  bins  than  the  900  s  measurement. 
The  pulse  amplitudes  did  not  exceed  400  mV,  and  thus  each  pulse  was  tallied  at  or 
below  channel  40.  A  constant  fluctuation  in  the  noise  baseline  was  observed  on  the 
oscilloscope,  most  likely  due  to  changes  in  the  leakage  current. 

4.1.1  UNL  resistive  device 

The  most  interesting  resistive  device  spectral  feature  was  the  broad  peak  centered 
at  approximately  channel  500.  In  Figure  28(a),  the  x-axis  corresponds  to  the  channel, 
or  energy  bin  number.  The  multichannel  analyzer  measures  each  incoming  pulse 
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Figure  28.  (a)  UNL  resistive  device  pulse  height  spectra  with  300/600/900  s  durations, 
(b)  The  number  of  total  counts  increased  and  the  dominant  peak  shifted  to  the  right 
100  channels,  meaning  the  amplitude  of  the  pulses  from  the  detector  had  increased 
with  cadmium  shielding. 


height,  assigning  the  height  to  one  of  many  pulse  height  ranges.  The  spectra  shown 
in  Figures  28(a)  and  28(b)  used  a  conversion  gain  setting  of  1024.  This  meant  the 
10.0  V  effective  range  of  the  multichannel  analyzer  was  divided  into  1024  equal  groups, 
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or  bins.  Once  the  incoming  pulse  was  properly  assigned  to  an  energy  bin.  a  counter 
incremented  the  corresponding  channel’s  count  number. 

The  peak  observed  was  first  believed  to  be  one,  or  two,  full  energy  peaks  from 
the  10B(n,  a)7Li  reactions.  Although  the  peak  was  somewhat  less  obvious  for  the 
300  s  duration  measurement,  it  was  plainly  visible  for  the  600  and  900  s  spectra.  The 
centroid  of  the  peak  shifted  slightly  for  the  longer  measurement,  but  its  shape  was 
well-preserved,  and  the  total  number  of  counts  under  the  peak  increased  for  the  longer 
measurements.  The  large  width  of  the  peak  was  consistent  with  lower  charge  carrier 
mobilities  expected  from  resistive  detectors  that  would  broaden  the  charge  collection 
time  distribution. 

The  device  was  wrapped  in  a  cadmium  foil,  and  a  pulse  height  spectrum  was 
collected  for  300  s.  The  counts  under  the  spectral  peak,  and  the  total  counts  for  the 
entire  spectrum  increased.  The  centroid  of  the  peak  also  shifted  to  a  higher  energy 
group.  These  three  observations  indicated  that  the  peaks  observed  were  induced  by 
something  other  than  a  neutron-induced  reaction. 

It  is  possible  that  thermal  neutron  interactions  occurred  in  spite  of  the  cadmium 
shielding.  Fast  neutrons  can  deposit  energy  in  the  foil  and  reach  thermal  equilibrium 
through  clastic  scattering  before  interacting  with  the  detector.  This  occurrence,  al¬ 
though  possible,  is  improbable  and  is  not  a  realistic  explanation  for  the  increase  in 
the  observed  count  rate.  This  increase  in  the  count  rates  and  peak  intensities  leads 
to  one  valid  conclusion.  Neutron  capture  reaction  products  are  not  responsible  for 
the  observed  pulse-height  spectrum  peak. 

Cadmium  absorbs  thermal  neutrons  by  (4)  and  emits  high  energy  y’s.  This  can 
cause  Compton  scattering  with  gadolinium  (acompton=22.7  b  [7]),  which  produces  low 
energy  gamma/X-rays  that  can  produce  pulses  within  the  resistive  bulk.  Another, 
and  potentially  more  likely  explanation  can  be  linked  the  the  PuBe  neutron  spec- 
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Figure  29.  Scope  trace  (infinite  persistence)  of  the  UNL  Resistive  Device  being  irradi¬ 
ated  by  PuBe  source. 

trum.  Elastic  and  inelastic  scattering  occurs  within  the  cadmium,  and  as  reported  by 
[44],  the  population  of  photons  below  20  keV  increases  due  to  down-scattering  in  the 
cadmium.  Because  lower  energy  photons  have  an  increased  probability  of  interaction 
within  the  detector,  this  is  a  logical  explanation  for  the  data  collected  from  this  re¬ 
sistive  device.  Additionally,  the  gadolinium  Compton  scatter  cross  section  for  20  keV 
photons  (o~compton=20.7  b  [7])  would  also  facilitate  an  influx  of  low-energy  photons 
capable  of  depositing  energy  within  the  detector. 

The  pulse  traces  were  observed  from  an  oscilloscope  (Figure  29)  during  the  irradi¬ 
ation  and  a  persistent  4.0  -  5.0  V  pulse  occurred  several  times  a  second.  The  device 
was  only  biased  to  90.0  V,  but  the  goal  was  to  apply  greater  voltages  (up  to  700  V)  in 
order  to  rapidly  sweep  any  neutron  induced  charges  towards  the  collecting  electrode. 
The  high  voltage  supply  was  set  to  the  maximum  output  yield  but  the  resistance 
across  the  preamplifier  circuit  greatly  attenuated  the  bias  strength.  The  apparent 
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noise  observed  on  the  oscilloscope  increased  with  bias  as  well.  Upon  considering  the 
pulse  height  spectra  collected,  improving  the  electrical  configuration  was  declined  due 
to  the  disadvantageous  effect  of  gamma/x-ray  interference. 

The  pulse  shown  is  characterized  by  two  features.  The  first  and  and  most  promi¬ 
nent  feature  is  the  large  amplitude  pulse,  with  a  short  rise  time.  This  is  characteristic 
of  the  fast  timing  properties  of  fully-depleted  semiconductor  detectors.  The  device 
was  resistive  and  depiction  was  not  expected.  These  large  peaks  with  short  rise 
times  are  likely  due  to  the  fluctuating  leakage  current  of  the  device.  The  pulse  tail 
shape  does  however  resemble  the  expected  pulse  shape.  The  high  resistivity  of  device 
would  naturally  reduce  the  charge  carrier  drift  velocity  and  increase  charge  collection 
times,  which  is  consistent  with  the  small-amplitude  but  broad  peak  observed  on  the 
oscilloscope. 

4.1.2  UMKC  diode 

The  UMKC  diode,  like  the  UNL  resistive  device,  experienced  some  interaction 
when  placed  over  the  PuBe  source.  Oscilloscope  pulse  traces  showed  a  700  mV  peak 
was  present  with  an  approximate  frequency  of  once  every  5  s.  Pulse  height  spectra 
on  the  other  hand  did  not  show  any  distinct  features,  and  thus  a  second  collection 
iteration  was  performed.  Pulse- height  spectra  for  600/1800/3600  s  were  collected  in 
order  to  increase  the  overall  number  of  counts,  and  the  conversion  gain  setting  on 
the  multichannel  buffer  was  reduced  to  512  channels.  The  linear  amplifier  gain  was 
also  adjusted  to  produce  a  5.0  V  output  pulse  from  the  preamplifier  (observed  on  the 
oscilloscope). 

Increasing  the  time  duration  for  the  spectra,  and  increasing  the  amplifier  gain  did 
not  result  in  improved  spectra  with  identifiable  features  (Figure  30(a)).  Cadmium 
shields  again  were  used  to  help  discriminate  between  interactions  that  were  neutron 
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Figure  30.  (a)  UMKC  diode  pulse  height  spectra  with  600/1800/3600  s  durations,  (b) 
No  change  observed  with  cadmium  foil. 

induced  and  those  that  were  not.  The  cadmium  foil  had  no  effect  on  the  spectrum 
(Figure  30(b)).  The  spectra  collected  with  and  without  the  cadmium  shielding  over¬ 
lapped,  meaning  the  total  number  of  counts  per  channel  remained  relatively  constant. 
The  absence  of  any  spectral  feature  however  could  not  alone  invalidate  the  hypothe- 
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sis  central  to  this  research.  Hetero-junction  diodes  containing  neutron-sensitive  con¬ 
stituents  may  be  able  to  detect  neutrons,  but  the  electronic  and  material  properties 
may  only  produce  signals  above  the  noise  threshold  with  precisely  optimized  geome¬ 
tries  and  fabrication  procedures. 

The  pulses  observed  in  the  oscilloscope  (Figure  31)  were  the  result  of  some  effect 
occurring  in  the  detector  induced  by  the  PuBe  source.  Background  oscilloscope  traces 
did  not  show  these  pulses  or  any  type  of  signal  other  than  small  fluctuations  in  the 
noise  baseline.  Persistence  was  enabled  (10  s  variable  persistence)  as  a  means  of 
comparing  the  pulse  amplitudes.  The  overlapping  peaks  shown  in  Figure  31  have 
amplitudes  ranging  from  2.0  -  5.0  V.  This  continuum  of  peak  heights  caused  higher 
channel  counts  on  the  pulse  height  spectra,  but  the  distribution  did  not  favor  any 
one  particular  energy.  Thus,  no  features  were  apparent  on  the  pulse  height  spectra 
in  spite  of  the  oscilloscope  output. 

The  leakage  current  increases  noise  baseline.  The  leakage  current  can  also  fluc¬ 
tuate,  and  this  is  known  source  of  noise  in  semiconductor  detectors  [14].  This  effect 
exists  in  both  bare  and  cadmium  shielded  spectra,  and  it  could  be  responsible  for  the 
observed  activity  on  the  oscilloscope. 

The  pulse  shape  of  partially-depleted  detectors  is  characterized  by  two  features. 
The  induced  depletion  current  contributes  to  the  fast  rise  time  and  large  pulse  am¬ 
plitude  [49].  Detectors  operated  at  or  above  saturation  have  constant  electric  fields 
throughout  the  device  and  carrier  transport  is  a  drift  process  only,  and  hence  the 
charge  collection  times  are  small.  In  partially  depleted  detectors,  the  charge  carrier 
transport  occurs  by  both  drift  and  diffusion  processes.  The  vanishing  electric  fields 
cause  a  reduced  drift  velocity  as  carriers  migrate  toward  the  electrode,  resulting  in 
drawn-out  pulse  tails  [16]. 

Recall  from  Figure  20(a),  the  UMKC  diode  was  depleted  less  than  10%  into  the 
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Figure  31.  Scope  trace  (infinite  persistence)  of  the  UMKC  diode  being  irradiated  by 
PuBe  source. 


silicon.  Pulse  tails  would  be  expected  with  this  limited  amount  of  depletion,  and  their 
absence  in  the  oscilloscope  output  was  a  clue  that  the  observed  pulses  on  the  scope 
were  not  related  to  ionized  charge  transport  through  the  detector. 

Neutron  capture  products  were  created  in  the  boron  carbide  region,  and  caused 
ionizations  in  the  depleted  silicon.  It  is  also  feasible  that  the  induced  signal  was 
below  the  minimum  detectable  limit  of  the  collection  electronics.  The  drift-induced 
contributions  to  the  current  pulse  were  adversely  affected  by  low  electric  fields.  Once 
charge  entered  the  un-depleted  bulk  region  of  the  silicon,  the  electric  held  was  no 
longer  present  and  carrier  transport  occurred  via  diffusion  only. 

The  maximum  charge  collection  length  in  partially  depleted  detectors  is  given  by 


Lp  — 


2Ldiff  + 


(26) 
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where  Ldiff  is  a  length  by  charge  diffusion  [50].  Large  diffusion  lengths  contribute 
to  longer  carrier  lifetimes  and  enhance  charge  collection  probabilities.  The  diffusion 
length  is  given  by  (29)  and  (27)  is  used  to  calculate  the  minority  carrier  mean  free 
time.  The  diffusion  length  in  n-type  silicon  is  determined  by  the  minority  charge 
carrier  type,  which  is  denoted  by  the  p  subscript  in  (27  -  29). 


Tp 

(27) 

(kT\ 

(28) 

Dp  —  1  —  J  pp 

Ldiff  =  \/Dptp 

(29) 

The  carrier  lifetime,  rp,  is  determined  from  the  hole  mobility  and  effective  mass, 
mp.  In  silicon,  the  effective  mass  0.69mo  where  mo=0.91x  10"30  kg.  The  diffusion 
coefficient,  Dp ,  is  determined  from  (28),  which  is  known  formally  as  the  Einstein 
relation  [15]. 

The  diffusion  length  calculation  (Appendix  B.1.1)  gives  a  diffusion  length  of 
1.4  /ini.  With  the  silicon  only  depleted  44.5  /im  from  the  junction,  (26)  yields  a  charge 
collection  length  slightly  less  than  25%  of  the  silicon  wafer  thickness  (Lp  «  180  pm). 
This  also  implies  the  induced  charge  decays  significantly  once  swept  through  the  de¬ 
pletion  region  and  probability  of  absorption  prior  to  being  collected  by  the  electrode 
is  great. 

Both  the  insufficient  depiction  width  and  total  charge  collection  length  (26)  were 
likely  causes  for  the  absence  of  any  detector  response.  A  simple  device  model  was 
created  to  analyze  the  detector  charge  collection  capabilities  and  efficiencies.  Geo¬ 
metric  parameters  of  the  diode,  such  as  the  silicon  and  boron  carbide  thicknesses, 
and  the  biases  applied  to  the  electrodes  were  varied  in  order  to  identify  factors  which 
contributed  to  the  negative  detector  response. 
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4.2  Diode  Modeling 


Davinci  is  a  device  simulation  program  that  models  the  behavior  of  wide  variety  of 
semiconductor  devices,  including  hetero-junction  diodes,  in  two  or  three  dimensions. 
The  software  numerically  solves  Poisson’s  equation  and  both  continuity  equations 
(electrons  and  holes)  in  order  to  analyze  the  device’s  electrical  characteristics  at  an 
arbitrary  bias.  Due  to  the  complexities  of  semiconductor  devices,  the  solutions  to 
this  system  of  three  partial  differential  equations  are  rarely  obtained  analytically. 
Analytical  solutions  are  also  only  achieved  after  several  simplifying  assumptions  have 
been  made.  Davinci  discretizes  the  three  equations  over  a  nodal  grid  and  produces 
accurate  solutions  to  them. 

The  specific  objective  of  the  modeling  was  to  replicate  the  irradiation  experi¬ 
ment  described  above.  The  depletion  width  and  saturation  current  were  affected  by 
varying  user-defined  parameters  in  order  identify  the  source  of  the  device’s  failure 
in  producing  a  signal.  Davinci  replicated  the  ionization  caused  by  energetic  charged 
particles  within  that  depletion  region,  and  produced  a  time-dependent  transient  so¬ 
lution.  Using  transient  solution  outputs,  the  induced  charge,  potential,  and  current 
were  calculated  at  the  collecting  electrodes.  The  initial  step  in  the  modeling  was  the 
creation  of  a  diode  and  mesh  that  replicated  the  I-V  characteristics  of  the  actual 
physical  diode.  This  was  then  followed  by  heavy  charged  particle  ionization  simula¬ 
tion.  The  depletion  widths  were  then  extended  further  into  the  silicon  by  increasing 
the  reverse  bias.  The  diode  geometric  parameters  were  also  varied  in  order  to  conduct 
a  small-scale  optimization  study  and  determine  which  factor  played  the  greatest  role 
in  causing  the  negative  detector  response. 
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4.2.1  Creating  the  diode  and  mesh 


Figure  32  illustrates  the  3D  representation  of  the  modeled  diode.  Only  a  small 
section  of  the  diode  could  be  replicated  in  Davinci  due  to  the  limited  number  of 
available  nodes.  The  width  of  the  diode  was  retained  in  order  to  give  fidelity  to 
model- generated  output.  The  suspected  issue  with  the  diodes  was  the  excessive  silicon 
width,  and  the  inability  to  adequately  deplete  it  or  deliver  charge  to  the  electrodes. 

Amorphous  boron  carbide  was  not  available  as  material  within  the  Davinci  li¬ 
brary.  The  generic  semiconductor  (SEMICOND)  was  used  as  the  base  material  and 
its  default  properties  were  modified  according  to  Table  6  so  that  the  diode  was  a 
realistic  representation  of  the  physical  diode  used  during  the  actual  experiment.  Na 
values  in  Table  6  were  based  on  [31]  conflict  with  Table  4,  whose  values  correspond 
to  calculations  performed  after  the  device  modeling  was  completed.  Despite  this 
discrepancy,  the  modeled  device  had  a  larger  depletion  region  than  was  determined 
from  experimental  data.  With  the  modeled  diode’s  enhanced  depletion  and  hence, 
greater  detector  volume,  the  simulation  of  heavy  charged  particle  ionization  never¬ 
theless  failed  to  represent  sufficient  charge  to  be  detected  above  intrinsic  noise  of  the 
Ortec  142IH  preamplifier.  For  this  reason,  the  conclusions  drawn  from  the  model  are 
still  considered  relevant  and  valid. 

4.2.2  Model  versus  physical  device  comparison 

The  simulated  diode  current  density  characteristics  are  shown  in  Figure  33(a). 
J-V  was  a  means  for  comparing  the  characteristics  of  two  geometrically  dissimilar 
diodes.  The  reverse  bias  current  density  in  the  model  approached  the  magnitude  of 
the  physical  diode  reverse  current  density  by  a  factor  of  2.3,  and  within  a  factor  of 
1.3  of  the  theoretically  computed  saturation  current  density  from  (8).  The  diode  also 
replicated  the  sub-linear  forward  current  behavior  discussed  in  Section  3.2.1.  At  low 
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Figure  32.  3D  plot  of  the  diode  created  in  Davinci.  The  500  x  500  x  500  fim 3  Davinci 
model  essentially  represents  just  a  small  section  of  the  physical  diode  (illustrated  by  a 
small  cube  extracted  from  the  center  of  the  aluminum  contact).  The  disc-like  surface 
represents  the  boron  carbide  deposition  on  the  n-type  silicon. 


Table  6.  Material  parameters  for  Davinci  model 


Parameter 

Silicon 

Boron  Carbide 

Semiconductor  Material 

Silicon 

SEMICOND 

Prohle 

N-Type  (1011) 

P-Type  (4.2 xlO12) 

Bandgap 

Default 

0.9  [eV] 

AiO 

10-7  [s] 

10-8  [s] 

TpO 

KT6  [s] 

10-8  [s] 

Permittivity  Constant 

Default 

8 

Affinity 

Default 

4.58  [36] 

Nc 

Default 

10 19  [cm’3] 

Nv 

Default 

10 19  [cm-3] 

biases,  the  modeled  diode  exhibited  immeasurably  small  currents.  However,  the  J-V 
characteristics  at  5  V  reverse  bias,  the  same  conditions  set  during  the  irradiation  phase 
of  this  experiment,  adequately  replicated  the  actual  diode.  Therefore,  it  was  inferred 
that  individual  parameters  of  this  model  could  be  adjusted,  and  corresponding  effects 
on  the  virtual  device  would  be  the  same,  at  least  qualitatively,  to  similar  modifications 
of  the  actual  diode  parameter. 

The  diode  electric  held  decreases  linearly  from  the  junction  (Figure  4).  The  held 
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Figure  33.  Davinci  simulation  of  (a)  current  density  (J— V)  and  (b)  electric  field  at  5 
V  reverse  bias.  The  inset  is  a  close-up  view  of  the  the  first  10  microns  of  the  diode 
(starting  from  the  boron  carbide  surface). 


extends  to  a  depth  that  contradicts  calculations  presented  in  Appendix  B.1.1  predict¬ 
ing  approximately  45  /irn  of  depletion  at  5  V  reverse  bias.  The  electric  field  is  initially 
zero  on  the  p-side  contact  and  rapidly  increases  to  the  maximum  field  strength  value 
of  2156  [V/cm]  at  the  junction.  From  the  peak  electric  field  at  the  junction  and 
extending  towards  the  silicon,  the  electric  held  becomes  discontinuous  and  drastically 
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decreases  in  strength  by  a  magnitude  of  1400  [V/cm]  within  the  first  4  microns  into 
the  silicon.  The  inset  in  Figure  33(b)  shows  this  discontinuity  caused  by  the  differ¬ 
ences  in  the  semiconductor  permittivities.  Following  the  sharp  decrease,  the  electric 
held  drops  linearly,  as  expected,  until  it  decays  almost  completely  at  300  /jrn  into  the 
silicon. 

4.2.3  Heavy  particle  ionization  simulation 

The  purpose  of  this  phase  of  the  modeling  was  to  observe  the  transient  output  on 
the  positive  (n-side)  electrode  once  a  virtual  heavy  charged  particle  was  injected  into 
the  silicon.  Some  simplifications  were  necessary  to  allow  for  the  simulation  to  proceed, 
and  without  excessive  run  times.  First,  the  particle  track  was  assumed  to  follow  a 
straight  line  parallel  to  the  y-axis  of  the  diode  and  originate  at  the  the  coordinates 
(x,  y,  z)  =  (0,  0,  0).  The  column  of  charged  particles  was  assumed  to  have  a  Gaussian 
radial  dependence.  The  temporal  generation  dependencies  were  also  assumed  to  be 
Gaussian. 

The  a  particle  had  the  greatest  kinetic  energy.  Therefore,  it  caused  the  most 
ionizations  and  with  the  longest  ionization  track  length  through  detector.  The  energy 
deposited  in  the  detector  was  thus  based  on  the  a  stopping  power  of  silicon.  Charged 
particle  injection  was  simulated  through  the  PHOTOGEN  command. 

User-defined  parameters  provide  the  ability  to  customize  the  simulation  and  con¬ 
trol  the  amount  of  charge  injected  and  from  where.  PHOTOGEN  also  provides  the 
option  to  call  a  linear  energy  transfer  hie  in  order  to  prescribe  the  amount  of  charge 
deposited  as  a  function  of  depth.  This  energy  transfer  table  originates  from  a  stop¬ 
ping  power  data  [51],  and  (30)  facilitates  a  conversion  to  pC//nn.  The  values  are  then 


67 


in  a  form  that  is  compatible  with  Davinci. 


LET 
Lf(l ) 


—  [MeV  cm2  g"1] 
ax 


rl  LET  p 
3  G  /? 


dx  [pC  /ini' 


-ii 


(30) 


In  (30),  p  is  the  density  of  silicon  (2.33  [g  cm"3]),  and  e,  is  the  ionization  energy  of 
silicon.  A  MATLAB  script  interpolates  the  table  values,  and  iteratively  decrements 
the  particle  energy  for  every  0.1  /mi  (dE  =  dE/dx  Ax)  so  that  the  linear  energy 
transfer  values  can  be  obtained  with  respect  the  depth  into  the  silicon.  The  output 
of  the  script  is  plotted  in  Figure  34. 


Figure  34.  Particle  energy  and  LET  in  silicon  as  a  function  of  depth.  The  left  ordinate 
corresponds  to  the  particle  energy  (shown  in  blue);  the  right  ordinate  corresponds  to 
the  LET  in  silicon. 


The  charge  generation  function  has  linear  dependency,  as  well  as  time  and  radial 
dependence.  Davinci  simulates  a  Gaussian  radial  and  temporal  charge  distribution 
by  (31),  and  (32)  respectively. 


R(r)  =  exp 


DCHR 


(31) 


In  (31),  r  is  the  radius  and  DC  HR  specifies  the  characteristic  radial  charge 
distance[52].  Simulation  using  the  TRIM  code  report  that  5  nm  is  the  character¬ 
istic  radial  charge  distance  of  1.47  MeV  a  particles  in  silicon  [36].  A  mesh  scheme 
with  such  a  fine  spacing  was  impossible  due  to  the  size  of  the  device  and  node  lim¬ 
itations  invoked  by  the  software.  The  characteristic  radial  distance  was  increased 
to  0.2  /jrn  in  order  to  prevent  excessively  long  run  times.  In  (32),  TO  and  TC  are 
the  offset  time  and  pulse  generation  time  respectively.  The  offset  time  is  arbitrarily 
chosen,  and  the  pulse  generation  time  is  1.9xl0"12  s  (calculated  from  stopping  power 
data). 


Figure  35  is  the  time-dependent  transient  solution  for  a  single  a  particle  depositing 
its  energy  within  the  active  region  of  the  detector.  The  current  shown  was  adjusted 
to  account  for  differences  between  the  virtual  and  the  physical  diode  cross-sectional 
areas.  The  magnitude  of  the  pulse  was  compared  to  the  estimated  leakage  current,  and 
the  pulse  was  also  integrated  to  calculate  the  total  charge  induced  on  the  electrodes. 
Both  of  these  calculations  were  simplified  by  correcting  the  current  to  account  for  the 
cross  sectional  area  differences  instead  of  using  the  current  density.  The  pulse  peaked 
in  less  than  10  ps,  before  decaying.  The  current  then  continued  increasing  until  it 
reached  the  steady-state  reverse  bias  current. 

Drift  velocity  is  proportional  to  the  electric  field  (15),  and  given  the  low  field 
strength  (Figure  33(b))  in  this  diode,  the  collection  time  is  orders  of  magnitude  larger 
than  initially  calculated.  The  transient  current  is  given  by  dQ/dt,  and  as  collection 
times  become  long,  the  transient  currents  become  smaller.  Another  point  of  consid¬ 
eration  is  the  enhanced  recombination  and  trapping  that  occurs  with  decreased  drift 
velocity. 

One  conclusive  measure  of  whether  or  not  the  pulse  was  detectable  was  to  com- 
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Figure  35.  Simulated  pulse  output  for  a  simulated  diode  using  the  Davinci  transient  so¬ 
lution  mode.  The  time-dependent  solution  is  computed  following  a  simulated  injection 
of  charge  via  the  PHOTOGEN  command. 


pare  the  total  induced  charge  with  the  intrinsic  noise  charge  of  the  charge-sensitive 
preamplifier  used  during  this  experiment.  The  preamplifier  intrinsic  noise  level, 

FWHM  (noise)  =  2.35e^^,  (33) 

Vp 

is  commonly  reported  as  a  function  of  input  capacitance.  For  the  case  of  this  exper¬ 
iment,  the  intrinsic  noise  level  was  4.6  keV[42],  e  is  the  amount  of  energy  required 
to  create  one  charge  carrier  (3.6  eV  for  silicon),  the  Vrms/Vp  is  the  intrinsic  noise 
charge[53].  The  intrinsic  noise  charge  calculated  via  (33)  reached  2.04xl0"16  C.  The 
total  induced  charge  was  calculated  by  (34)  where  Q  was  the  total  induced  charge,  tc 
was  the  charge  collection  time,  and  current  I(t)  was  a  function  of  time. 

Q=  I  I(t)dt  (34) 

Jo 

The  total  induced  charge  (4.67xl0-1'  C)  did  not  exceed  the  intrinsic  noise  charge  of 
the  preamplilihcr.  This  meant  the  charge  detection  capability  of  the  preamplifier  was 
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(a) 


(b) 


Figure  36.  Davinci  simulation  showing  the  effect  of  bias  on  (a)  the  electric  field  (the 
inset  figure  is  an  expansion  of  the  top  3  microns)  and  (b)  the  transient  current. 

not  sensitive  to  the  quantity  of  induced  charge  resulting  from  the  neutron  capture. 
Figures  36(a)  and  36(b)  investigate  the  device’s  electric  fields  and  pulse  outputs  with 
100  microns  of  silicon  removed.  Additionally,  the  width  of  the  boron  carbide  was 
doubled.  Increasing  the  amount  of  boron  carbide  theoretically  would  increase  the 
degree  of  depletion  in  the  silicon  region,  but  the  opposite  effect  resulted  from  the 
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geometry  just  described.  The  electric  fields  decreased  as  a  result  of  increasing  the 
boron  carbide  thickness.  The  pulses  generated  from  the  Davinci  transient  solution 
showed  a  suppressed  peak  at  5  V  reverse  bias.  The  pulse  did  not  exceed  the  steady 
state  saturation  current  with  10  V  reverse  bias  applied.  Boron  carbide  and  silicon 


Depth  [|im] 

(a) 


(b) 


Figure  37.  (a)  Electric  held  and  (b)  transient  current  with  only  10  microns  of  silicon. 


effective  doping  concentrations,  the  particle  strike  origin,  and  particle  direction  were 


also  examined  to  determine  what  effect,  if  any,  they  had  on  the  detector  performance. 
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After  a  completing  multiple  simulations,  the  electric  field  strength  was  the  only  factor 
that  noticeably  affected  the  diode  charge  collection  ability  and  the  resultant  transient 
current.  This  conclusion  was  supported  in  Figures  37(a)  and  37(b).  The  diode  used 
in  this  simulation  contained  a  1  micron  thick  boron  carbide  region  and  a  10  micron 
silicon  region. 

This  diode  was  depleted  beyond  the  saturation  point,  and  the  electric  held  was 
large  and  constant  throughout  the  silicon.  The  transient  solutions  were  normalized 
in  order  to  account  for  the  different  steady  state  currents  at  5.0  -  20.0  V  reverse 
bias.  At  5.0  V  reverse  bias,  the  transient  current  pulse  was  1.78  times  greater  than 
the  reverse  bias  current;  at  20.0  V  reverse  bias,  the  transient  maxima  reached  2.25 
times  the  magnitude  of  the  steady-state  current,  ft  was  clear  that  electric  fields  in  the 
UMKC  diode  were  too  weak  and  the  resultant  slow  drift  velocities  hindered  the  charge 
collection  efficiency  by  increasing  the  charge  collection  times  and  the  probability  of 
recombination  and  trapping. 
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V.  Conclusion 


5.1  Experimental  Findings 

No  experimental  evidence  of  neutron  detection  was  observed  during  the  course 
of  this  experiment.  The  devices  were  introduced  into  thermal  neutron  environments 
with  a  calculated  total  neutron  flux  of  464  [n/cm2-s].  Neutrons  interacted  within  the 
diodes,  but  no  corresponding  signal  was  identified  on  either  the  oscilloscope  or  the 
pulse-height  spectra.  The  root  cause  of  a  failed  detector  response  was  investigated 
through  modeling. 

The  following  summarizes  the  findings  from  this  research: 

1.  The  UNL  resistive  device  exhibited  behavior  entirely  opposite  to  the  anticipated 
response  to  shielding  the  device  with  a  cadmium  foil.  The  pulse  magnitudes  and 
intensities  increased  with  the  addition  of  the  cadmium  shield.  This  was  a  clear  indi¬ 
cation  that  the  observed  peak  in  the  pulse-height  spectra  was  caused  by  interactions 
other  than  thermal  neutron  capture.  The  cadmium-covered  device  also  demonstrated 
the  device’s  gamma  sensitivity. 

Cadmium  is  a  high  Z  element  and  can  soften  the  PuBe  photon  spectrum  by  down¬ 
scattering  and  the  result  is  a  higher  flux  of  photons  below  20  keV.  Since  lower  energy 
photons  have  a  higher  probability  of  interacting  with  the  detector,  it  is  likely  that  this 
was  a  source  of  photon  interference  observed  in  the  oscilloscope  while  the  detector 
was  being  irradiated.  Additionally,  cadmium  absorbs  thermal  neutrons  and  re-cmits 
263.5  keV  gamma  rays  as  a  product  of  the  radiative  capture  reaction.  Gadolinium, 
which  is  also  a  high  Z  material,  has  a  Compton  scatter  cross  section  of  20.1  b  and 
22.7  b  for  20.0  and  263.5  keV  photons  respectively.  The  gadolinium  oxide  film  may 
have  enhanced  the  device’s  photon  sensitivity  by  down-scattering  the  photons  and 
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increasing  the  population  of  photons  depositing  energy  within  the  detector. 

2.  The  UNL  diode  junction  was  likely  damaged  during  the  course  of  experiments 
performed  previously.  Useable  pulse-height  spectra  could  not  be  obtained  from  the 
diode.  Pre-characterization  of  the  device  showed  the  device’s  C-V  characteristics 
deviated  from  the  expected  junction  behavior.  This  was  the  first  indication  that 
problems  persisted  within  the  diode,  and  it  was  not  surprising  that  diode  performed 
poorly  while  being  irradiated.  The  dopant  concentrations  N a  and  Np  could  not  be 
calculated  from  the  measured  C-V  data,  and  thus  a  theoretical  reverse  saturation  cur¬ 
rent  density  was  not  approximated.  The  diode  was  less  conductive  than  the  UMKC 
diode,  but  without  an  estimate  for  the  reverse  saturation  current  density,  the  leakage 
current  remained  unknown,  but  appeared  to  dominate. 

3.  The  UMKC  diode  exhibited  large  leakage  currents  (~10  fiA  estimated)  at  5  V 
reverse  bias.  This  exceeded  the  ideal  by  4  orders  of  magnitude.  C-V  measurements 
indicated  that  the  silicon  was  being  depleted,  but  the  depletion  width  did  not  extend 
beyond  25%  of  the  silicon  depth.  The  built-in  voltage  and  dopant  concentrations 
were  calculated  with  the  pre-characterization  data,  but  the  silicon  dopant  concentra¬ 
tion  computed  was  two  orders  of  magnitude  higher  than  expected.  This  may  have 
contributed  to  high  capacitance,  which  adversely  affected  the  S/N  ratio,  and  caused 
electric  fields  within  the  diode  to  be  lower  than  ideal. 

Recombination  rates  in  the  depletion  region  are  improved  with  slow  ion  drift  ve¬ 
locities.  The  weak  electric  fields  were  expected  given  the  depletion  width  calculated 
from  the  C-V  data.  Pre-irradiation  theoretical  calculations  suggested  that  the  tran¬ 
sient  current  collected  from  the  n-type  silicon  side  of  the  diode  would  exceed  the 
measured  steady  state  current,  but  only  by  a  small  margin. 
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The  detector  pulse-height  spectrum  presented  no  useable  features,  and  although 
the  oscilloscope  trace  activity  increased  in  activity  when  the  device  was  irradiated 
with  a  PuBe  source,  there  was  no  change  in  the  signal  between  the  unshielded  and 
cadmium-shielded  spectra.  Therefore,  nothing  observed  in  the  oscilloscope  could  be 
attributed  to  ionization  produced  from  the  10B(n,  a) 'Li  reaction.  Determining  the 
root  cause  for  this  in  the  single  pulse-height  spectrum  was  the  primary  objective  of 
device  modeling.  Experimental  conclusions  asserted  that  the  detector  response  would 
not  exceed  the  noise  background  given  the  weak  electric  fields  in  the  depletion  region. 
This  argument  was  supported  through  modeling. 

The  UMKC  diode  was  replicated  in  Davinci,  and  like  the  physical  diode,  simu¬ 
lation  produced  weak  electric  fields  in  the  depletion  region.  This  caused  the  charge 
collection  times  to  become  excessively  long  with  respect  to  the  fast  timing  character¬ 
istics  expected  of  semiconductor  detectors.  The  weak  electric  fields  also  led  to  higher 
recombination  rates  which  depleted  the  charge  carrier  population.  Ultimately,  it  was 
determined  that  the  charge  generated  by  a  single  a  particle  ionization  track  did  not 
exceed  the  intrinsic  noise  charge  of  the  preamplifier  used  for  this  experiment. 

The  effect  of  variations  on  other  parameters,  such  as  the  silicon  and  boron  carbide 
thickness,  carrier  concentrations,  and  ionizing  particle  paths,  were  analyzed.  When 
the  silicon  thickness  was  reduced  to  10  /im,  the  diode  was  over-depleted  with  only  5  V 
reverse  bias  and  the  electric  field  strength  was  constant  throughout  the  diode.  The 
charge  induced  with  this  geometry  produced  transient  current  pulses  1.75  times  the 
reverse  bias  steady-state  current.  The  conclusion  drawn  from  the  simulation  study 
was  that  thinner  silicon  regions  resulted  in  lower  saturation  voltages,  which  led  to 
higher  electric  fields,  improved  charge  carrier  recombination  escape  probabilities,  and 
thus  enhanced  the  transient  currents. 
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5.2  Recommendations  for  future  work 


The  device  modeling  supported  the  conclusion  that  charge  carriers  were  lost  due 
to  slow  drift  velocities,  but  this  complication  was  resolved  by  reducing  the  silicon 
thickness  so  that  the  diode  could  be  fully  depleted.  The  motivation  for  this  neutron 
detection  study  was  to  develop  light-weight  but  efficient  thermal  neutron  detectors 
that  could  be  used  in  ports  of  entry  to  detect  SNM  smuggling  attempts.  Diodes  such 
as  the  one  modeled  with  only  10  /im  would  be  worth  investigating  due  to  their  inher¬ 
ently  low  saturation  voltage.  These  diodes  would  likely  require  different  fabrication 
techniques  however  due  to  the  fragility  of  extremely  thin  silicon  wafers. 

Another  option  to  explore  would  be  the  use  of  silicon  wafers  half  the  thickness  of 
the  ones  used  during  this  study.  The  simulation  showed  the  depletion  region  depth 
extending  past  250  /an  at  5  V  reverse  bias.  The  detector  would  not  then  be  operated 
in  an  over-depleted  state,  but  it  could  possibly  still  collect  an  amount  of  charge  that 
would  exceed  the  intrinsic  noise  of  the  preamplifier.  Recall  also  that  the  intrinsic 
noise  level  is  often  a  function  of  capacitance.  The  capacitances  of  the  UMKC  diode 
was  large,  and  the  calculated  silicon  dopant  concentration  was  greater  than  expected. 
Improving  the  fabrication  techniques,  even  while  still  using  silicon  and  boron  carbide 
thicknesses  identical  to  the  diodes  used  in  this  study,  may  help  reduce  the  diode 
capacitances  and  therefore  the  charge-sensitive  preamplifier  intrinsic  noise  charge. 

The  PECVD  reactor  at  UMKC  was  recently  disassembled  and  rebuilt  to  facilitate 
an  improved  and  completely  in  vacuo  diode  fabrication  process.  The  silicon  wafers 
would  be  loaded  into  the  the  multi-chamber  reactor,  boron  carbide  and  contacts  would 
then  be  deposited  (or  sputtered)  without  exposing  the  devices  to  oxygen  or  other 
surface  contaminants.  Diodes  produced  via  this  procedure  will  likely  have  improved 
junctions,  both  at  the  p-n  junction  and  the  metal/semiconductor  junctions  where  the 
contacts  are  sputtered.  As  material  science  technologies  advance,  it  is  highly  likely 
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that  true  direct-conversion  detectors  could  be  produced.  It  is  well  worth  the  effort  to 
test  and  analyze  the  improved  diode  electronic  characteristics  and  any  pulse-height 
spectra  that  could  be  obtained. 

This  experiment  suffered  from  a  dearth  of  useable  diodes.  It  would  be  enlight¬ 
ening  to  acquire  a  set  of  samples,  all  created  during  the  same  deposition  process, 
and  compare  the  C-V  characteristics.  This  would  enable  a  comparison  between  the 
calculated  dopant  concentrations  in  order  to  determine  whether  the  high  dopant  con¬ 
centration  calculated  in  this  study  was  anomalous,  or  an  indication  of  some  other 
phenomena.  The  effective  dopant  concentration  is  of  fundamental  importance  and  it 
affects  the  depletion  widths  that  could  be  obtained.  It  is  thus  a  non-trivial  matter 
and  well-worth  investigation. 

Research  focused  on  the  trap  densities  and  boron  carbide/silicon  interface  defect 
densities  and  the  corresponding  effects  on  detector  efficiencies  may  result  in  method¬ 
ologies  that  profoundly  improve  the  performance  of  these  devices.  There  is  little 
agreement  between  the  electrical  property  values  of  boron  carbide  as  reported  by  dif¬ 
ferent  literature  sources.  The  modeling  effort  in  this  study  could  easily  be  augmented 
and  improved  with  corrected  boron  carbide  parameters.  Many  questions  were  left 
unanswered  by  the  simulations  including  the  topic  of  the  charge  collection  length  in 
partially-depleted  detectors,  and  the  effects  of  trapping.  Resolving  these  matters  may 
help  in  adding  fidelity  and  address  the  literature  differences  concerning  boron  carbide 
and  its  electrical  properties. 
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Appendix  A.  Equipment  and  Settings 


List  of  Electronic  Components 

Component 

Manufacturer 

Model 

Serial 

Pulser 

Ortec 

480 

2391 

Preamplifier 

Ortec 

142IH 

2035 

Linear  Amplifier 

Ortec 

572A 

536 

ADCAM  (MCB) 

Ortec 

926 

06165185 

Oscilloscope 

Tektronics 

DP07104 

B022189 

Bias  Supply 

Ortec 

478 

41975 

Summary  of  Component  Settings 

Component 

Parameter 

Diode  Device 

Resistive  Device 

Linear  Amplifier 

Shaping  Time 
Coarse  Gain 
Fine  Gain 

2  /xs 

500 

5.0 

2  /xs 

50 

5.65 

ADCAM  (MCB) 

LLD 

Conversion  Gain 

0.576  V 

512 

0.08  V 

1024 

Bias  Supply 

40  V 

1000  V 
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Appendix  B.  Calculations 


B.l  Diode  Parameter  Calculation 
B.1.1  UMKC  Built-in  Voltage 


1/C2,  given  the  linear  dependence  of  (Vh  -  V),  was  used  to  calculate  Vbi. 
Line  of  Fit:  y  =  6.929  x  1017  -  2.905  x  1018 V 


f±Y  =  f  2  ( eslND  +  ebcNA)  (Vhi  -  V) 
\Cj  V  qesiebcNDNA 

when  Vbi  =  V 


c)  ° 

0  —  m(V  —  Vbi )  +  b 
V^  =  0.239  [V] 


Nd  and  NA  could  then  be  determined.  Mathematica  “FindFit”  command  was  used 
to  accomplish  this. 


3.13  xlO13 

cm  3 

NA(bc ) 

7.0  xlO12 

cm  3] 

The  constants  used  for  the  calculations  above  are  given  below. 

q  =  1.602  xl0~19 
m  =  -2.905  xlO18 
b  =  6.929  xlO17 
k  =  11.9 

e0  =  8.85418  xlO’14 
es  =  ke0 
£bc  =  8e0 

The  UNL  diode  did  not  have  1/C2  linearity;  the  dopant  concentrations  and  built 
in  voltage  could  not  be  calculated,  but  the  capacitance  data  could  be  used  to  approx¬ 
imate  the  depletion  width  (shown  below). 
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B.1.2  UMKC  and  UNL  diode  hetero-junction  calculations 


Cj 


W  =  en 


+  xp 


UMKC 

UNL 

xp  [n m 

1.5 

2.0 

Depletion  width  approximation 

W(V  =  0) 

10.94 

136.1 

W(V  —  5) 

44.53 

128.8 

B.2  Pulse  Prediction  Calculations 


The  constants  used  to  calculate  charge  collection  times  (tc),  transient  currents  (Itransient) , 
and  time  constant  (r)  are  given  below. 

Edev  =  1.77xl06  [eV] 

<,=  1.602  x  10-i»  [C] 

Epr  =  3.62  [eV] 

t  plasma  5x10  [s] 

W  =  44.5261  x  10-4  [cm]  (UMKC  diode  depletion  width  at  5  V) 

/i Si  =  1400  [cm2/V-s] 

I^bc  =  0.00075  [cm2/V-s] 


Charge  collection  time 


tc 

S(x) 

tc(x) 


_d_ 

£fi 

qND 


(W-x) 


log 


W 


qfiND  W  —  Xq 
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tr  Si  [si  — > 


2.83  xlO-9 

s 

tc  BC  [s]— > 

1.78  xlO"4 

s] 

Charge  collection  times  can  differ  by  5  orders  of  magnitude.  Only  the  charge  collected 
within  the  instrumental  circuit  time  constant  will  be  integrated  on  the  preamplifier 
capacitor. 


Time  constant 


r 

C 


RC 

r  i 


R  —  93  [Q]  (reported  by  [42]) 

Cl  =  10000  xlO-12  [F]  (reported  by  [42]) 
C2  =  0.05 xlO"9  [F] 

Time  Constant,  r  — > 


measured) 


4.63x10 


-9 


Transient  current  calculation 


*■  transient 


=  2 


RdepQ 


Epr  (tc  T  t plasma) 


I  trails  (Si)  ^  2.00  XlO-5  [A]  Itrans(BC)— )■  8.79  X  10"10  [A] 


The  transient  current  for  the  Silicon  side  is  on  the  same  order  of  magnitude  as  the 
leakage  current  (estimated  10  /j,A)  at  the  5  V  reverse  bias.  The  transient  current  on 
the  BC  side  is  5  orders  of  magnitude  less  than  the  leakage  current.  Only  current 
collected  from  the  n-type  silicon  side  has  chance  of  being  observed. 


Diffusion  Length  calculation 

y/ '  DpTp 
mn 

q 


DP  =  (f )  Up 
k  =  1.38 xlO"23  [J/K] 
T  =  300  [K] 


Carrier  Lifetime  — > 


1.8xl0"9  [si  Diffusion  Constant 


11.6 


cm 


!/s 
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Diffusion  Length  — *  1.4  [fi m 


Saturation  Current 

q  =  1.602  xl0“19  [C] 

Nc  =  10 19  [cm-3] 

Ny  =  1019  [cm-3] 

Eg  =  0.9  [eV] 
k  =  8.6174  xl0~5  [eV/K] 

T  =  300  [K] 

Vth  =  f  [V] 

Silicon 

fie  =  1350  [cm2/V-s] 

Hh  =  450  [cm2/V-s] 

Boron  Carbide 
/ie  =  .001  [cm2/V-s] 

Hh  =  .00075  [cm2/V-s] 

Na  and  NB  were  calculated  from  C-V  data. 

Na  =  7.0  x  1012  [cm-3]  —  effective  dopant  concentration  in  boron  carbide 

Nd  =  3.13  x  1012  [cm-3]  —  impurity  concentration  in  silicon 

nil  =  0.69  mo  [kg]  —  mass  of  holes  in  silicon 

m2  =  5  mo  [kg]  —  mass  of  electrons  in  boron  carbide 

mo  =  0.91  xlO-30  [kg]  —  mass  of  electron 
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Js  =  Jn (Boron  Carbide)  +  Jp(Silicon) 
=  Dn'njq  Dprfq 
LuNa  LpNd 

=  5.3  xlO-5  [A/cm2] 


B.3  Neutron  Emission  Rate 


This  calculation  was  originally  performed  by  [45]  and  has  been  modified  according  to 
[43]  to  accurately  predict  to  the  neutron  emission  rate  at  the  time  of  the  experiment. 
NAv=  6.022  xlO23 

ti=24110 

t2=6564 

t3=14.35 

t4=432.2 


MassTot=76.36 
Wi=.9102 
W2=0. 08098 
W3=0. 00844 

Massi=MassTot  W4 
Mass2=MassTot  W2 
Mass3=MassTot  W3 

Atomic  Masses 
AM4=239. 052156 
AM2=240. 053808 
AM3=241. 0568453 


iVi=l;3 

\=\A 


Mass* 
AM  i 
Log  [2] 

2l 


NAv 


DateDifference[{1961,3,23},{1965,5,17},  “Year”] 
Time  elapsed  =  4.15068,  Years 
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This  time  will  be  used  to  calculate  Aq  within  the  exponential  terms. 


A 

Aq 

A[t] 


9.86  xlO6 


Ai 


1  +  E23  ((A.Nl  +  A^Oa-X) <Exp  1-4'151AJ  -  Exp  l-415A»]))) 

A4A3N3 


An  1  +  1.23 


(AjNl  +  A2N2)  (A3  -  A4) 


(Exp  [— A4t]  —  Exp  [— A3t]) 


v4[0]  =9.31466  x  106 

DateDifference[{1961,3,23},{2011,12,8},  “Year”] 
Time  elapsed  =  50.7104,  Years 


,4[50.7]  -)• 


1. 19182 xlO7  Is"1 


B.3.1  Reaction  Rate  Density 

^l.ig^xio7^-1] 

p=2.55  [g/cm3]  density  of  boron  carbide  [36] 
ma=78.161  [g/mol]  atomic  mass  of  boron  carbide  (B5CH12) 

Nav— 6.022x1023  [atoms/mol] 

<7^=3880  x  10"24  [cm2] 

r=0.252  [cm]  diode  radius  (assuming  spherical  geometry  due  to  isotropic  source  emis¬ 
sion) 

<7=20.32  [cm]  distance  between  detector  and  source 
B=. 94  branching  ratio  of  the  reaction 


RRD  = 

$th  = 


ABVl 

An 


D 

E 

N 


2n 


1  - 


\/r 2  +  d 2 
d 


CthN 

pmaNAV 


Thermal  Flux,  <5>th  — >  464  [n/cm2s 


Macroscopic  Cross  Section,  E  — > 

24578 

cm 

1 

Rate  Reaction  Density,  RRD  — > 

1.14xl07 

_ii 

cm  s  J 

85 


Appendix  C.  Davinci  Input 


TITLE  Diode  2  Simple  Diode 


MESH  DIAG.FL 

X.MESH  X.MIN=0  X.MAX=30.0  Hl=0.2  H2=10.  H3=0.2  SUMMARY 

X.MESH  X.MIN=30  X.MAX=250.0  Hl=10.  H2=100  H3=10.  SUMMARY 


Y.MESH  Y.MIN=0.  Y.MAX=10.0  Hl=0.1  H2=.5  H3=0.1  SUMMARY 

Y. MESH  Y.MIN=10.0  Y.MAX=500  Hl=0.5  H2=100  H3=.5  SUMMARY 

Z. MESH  Z.MIN=0  Z.MAX=30  Hl=0.2  H2=10  H3=0.2  SUMMARY 

Z.MESH  Z.MIN=30  Z.MAX=250  Hl=10  H2=100  H3=10  SUMMARY 


SDefining  the  region  materials  and  electrodes 

REGION  NAME=NREG  X.MAX=250.  Y.MIN=1.0  Y.MAX=500  Z.MAX=250.  SIL¬ 
ICON 

REGION  NAME=PREG  X.MAX=250  Y.MAX=1.0  Z.MAX=250  SEMICOND 
ELECTR  NAME=PCON  TOP 
ELECTR  NAME=NCON  BOTTOM 

PROFILE  REGION =NREG  UNIFORM  N-TYPE  N.PEAK=lell 
PROFILE  REGION =PREG  UNIFORM  P-TYPE  N.PEAK=4.2el3 

MATERIAL  TAUN0=lE-7  TAUP0=lE-6  SILICON 
SMOBILITY  SILICON  MUN0=0.001 
SMOBILITY  SILICON  MUP0=0.00075 

MATERIAL  SEMICOND  EG.MODEL=3  EG300=0.9  TAUN0=le-8  TAUP0=le-8 
PERMITTI=8 

+  AFFINITY=4.58  NC300=1E19  NV300=1E19 
SAVE  MESHFILE  OUT.FILE=DIODE2B.tdr  tdr 
MODELS  SRH  AUGER  FLDMOB  BGN  CCSMOB 

SCONTACT  NAME=NCON  RESIST=4.0E3  CAPAC=0.25E-13  PRINT 

SYMBOL  NEWTON  CARRIERS=2 
SOLVE  INIT 
METHOD  ITLIMIT=10 

PLOT. 3D  BOX  GRID  FILL  TITLE=  “PHOTODIODE  GRID  REGRID”  SCALE 
GAMER  A= (-500,-500,-500) 

SOLVE  V(NCON)=0  ELECTRODE=N CON  VSTEP=1  NSTEP=5 


PLOT. ID  E. FIELD  X.START=0  X.END=0  Y.START=0  Y.END=500  Z.START=0 
Z.END=0  TITLE=“ELEC  FIEL”  OUT.FILE=EFLD5V.dat 

PHOTOGEN  DCHR  =0.2  T0=3.0E-12  TC=1.9E-12  LETFILE=dvel821et  PC.UNITS 
$+  X.START=250.0  X.END=250.0  Y.MIN=1.  Y.MAX=10. 

$+  Z.START=250.  Z.END=250. 

SOLVE  TSTEP=0.5E-12  TSTOP=1.0E-10  OUT.FILE=alpha5Va 
PLOT. ID  X.AXIS=TIME  Y.AXIS=I(NCON)  BOTTOM=-2E-3 

+  COLOR=2  POINTS  TITLE=  “Example  4  -  Current  vs.  Time”  OUT.File=diode2BNPULSE.dat 
EXIT 
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